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ABSTRACT 


An  approach  to  solving  the  problems  encountered  in  dual 
quality  control  of  distillation  columns  is  tested  experi¬ 
mentally  using  a  pilot  scale  distillation  column  with  a 
methanol-water  system.  The  problems  dealt  with  include 
interaction  between  the  overhead  and  bottoms  composition 
control  loops,  measurement  noise  and  sluggish  response  of 
the  controlled  variables. 

To  lessen  these  problems  first  the  physical  configura¬ 
tion  was  studied.  An  experimental  program  was  carried  out 
to  evaluate  which  control  variable  for  overhead  composition 
control, ref lux  or  distillate  flow, gave  better  performance. 

It  was  found  that  depending  on  the  type  of  disturbance 
either  could  be  better  for  the  case  where  bottoms  composition 
control  was  not  operative.  The  distillate  flow  was  chosen  as 
the  control  variable  for  the  subsequent  control  studies  where 
both  terminal  compositions  were  regulated. 

A  control  technique  that  requires  relatively  simple 
modelling  to  obtain  decoupling  to  reduce  interaction  and 
feedforward  control  to  speed  disturbance  attenuation  was 
chosen  for  this  study.  A  linear  model  is  required  for  use 
in  the  formulation  of  the  control  law. 

Distillation  modelling  was  examined  with  respect  to 
mechanistic  type  models  using  energy  and  mass  balances  and 
physical  data.  A  concentration  model  was  chosen  for  the 
control  calculations  as  it  was  found  to  agree  with  experi- 
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mental  results  from  the  equipment  better  than  an  enthalpy 
model  that  was  tried. 

The  final  control  law  is  calculated  using  a  model 
reduced  from  the  original  tenth  order  concentration  model  to 
fourth  order  by  modal  reduction.  This  was  discretized  to 
enable  implementation  with  a  digital  computer. 

The  performance  of  the  system  was  simulated  and  the 
results  predict  good  disturbance  attenuation. 

Experimental  results  are  presented  for  the  control  tests 
with  the  equipment  including  those  in  which  integral  feedback 
was  added.  The  decoupling  in  the  control  function  facili¬ 
tates  on-line  tuning  of  feedback  functions.  Test  results 
presented  indicate  that  the  control  scheme  is  successful  for 
regulatory  control. 
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CHAPTER  1 


INTRODUCTION 


Control  of  a  binary  distillation  column  entails 
several  problems  that  make  simple  single  variable  feedback 
control  unsatisfactory.  First,  the  terminal  compositions 
respond  slowly  to  disturbances  in  the  feed  which  means  that 
a  large  upset  is  present  in  the  internal  variables  before 
control  action  can  be  taken.  This  gives  slow  recovery  from 
a  disturbance.  Second,  the  system  exhibits  highly  non-linear 
behavior  so  that  control  gains  must  be  changed  to  allow  for 
the  change  in  process  gain  when  the  system  is  away  from  the 
study  state  due  to  an  upset.  This  also  means  that  a  linear 
model  is  not  adequate  to  predict  the  response  of  the  system 
for  a  large  change  in  input  conditions.  Another  difficulty 
with  single  variable  feedback  control  stems  from  noise  in 
the  system  or  unreliable  measurements  which  may  result  in 
incorrect  control  action.  The  final  unfavorable  phenomenon 
is  the  interaction  that  takes  place  between  the  composition 
control  loops  on  the  overhead  and  bottoms  streams.  This 
reduces  the  stability  of  the  control  loops ,  degrades  the 
performance  and  makes  it  extremely  difficult  to  tune  the 
controllers  as  tuning  of  one  affects  the  other. 

The  object  of  this  work  is  to  develop  a  practical 
control  scheme  for  a  binary  distillation  column  for  regula¬ 
tion  of  the  terminal  compositions  and  to  present  experi- 
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mental  results  using  the  scheme  on  the  pilot  scale  distil¬ 
lation  column  at  the  University  of  Alberta.  The  development 
is  carried  out  in  three  stages,  namely:  (1)  an  evaluation 
of  the  choice  of  manipulated  variable  for  control  of  the 
overhead  composition,  (2)  modelling  of  the  distillation 
column,  and  (3)  the  calculation  and  implementation  of  the 
control  law. 

The  first  stage  of  this  work  is  an  experimental  study 
in  which  the  reflux  flow  and  the  distillate  flow  were  tested 
as  the  manipulated  variable  in  order  to  determine  which 
gives  better  control  of  overhead  composition.  Past  control 
studies  at  the  University  of  Alberta  on  the  column  have  used 
the  reflux  flow  exclusively,  and,  as  this  method  has  come 
into  question  (5,19,27,33)  it  is  compared  experimentally  to 
the  alternative  distillate  flow  for  the  case  where  only  the 
overhead  composition  is  controlled  while  the  bottoms  compo¬ 
sition  controller  is  left  on  manual. 

The  second  phase  of  this  study  is  the  derivation  of  a 
model  suitable  for  use  in  establishing  a  control  law.  The 
control  method  developed  in  this  study  is  restricted  to  a 
linear  model.  A  model  in  linear  state  space  format  which 
meets  that  criterion  was  developed  by  McGinnis  (28).  He 
used  a  twentieth  order  model  with  the  states  being  the 
liquid  enthalpies  and  flowrates  based  on  the  dynamic  energy 
balance  and  assuming  a  first  order  lag  for  the  liquid  flow- 
rates  from  tray  to  tray.  He  assumed  that  the  vapor  flow 
remained  constant  and  that  the  condenser  load  was  a  free 
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input.  This  model  gave  reasonable  results  when  the  experi¬ 
mental  change  in  the  condenser  load  was  added  to  the  dis¬ 
turbance  input  to  the  simulation.  It  should  be  noted  that 
the  condenser  load  is  used  to  control  the  column  pressure 
which  is  assumed  constant  so  the  condenser  load  is  deter¬ 
mined  by  the  other  inputs.  For  this  reason  other  methods  of 
representing  the  vapor  flow  and  condenser  load  were  con¬ 
sidered.  Non-linear  models  based  on  the  unsteady  state 
material  and  energy  balances  are  examined  along  with  the 
assumptions  necessary  to  make  them  solvable.  The  simulation 
results  from  the  models  are  then  compared  to  experimental 
results.  The  model  that  best  predicted  experimental  results 
was  linearized  around  steady  state  conditions  in  order  to 
enable  its  use  in  the  control  of  the  column. 

The  final  portion  of  the  study  is  devoted  to  obtaining 
a  control  scheme  that  can  be  implemented  on  the  column  and 
to  testing  it  experimentally.  The  control  law  was  based  on 
the  method  proposed  by  Greenfield  and  Ward  (13,14).  This 
method  gives  feedforward  control  with  input  and  output  de¬ 
coupling.  The  decoupling  removes  the  interaction  and  allows 
the  addition  of  simple  feedback  control  loops.  These  loops 
can  be  easily  tuned  on-line  so  that  the  errors  from  model¬ 
ling  and  linearization  can  be  compensated  for.  Berry  (2) 
developed  a  decoupled  feedback  control  scheme  that  gave 
excellent  control.  He  started  with  an  experimental  first 
order  plus  time  delay  transfer  function  model  but  found  that 
tuning  was  required  for  both  feedback  and  decoupling 
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elements.  McGinnis’  multivariable  scheme,  which  was  based 
on  a  second  order  discrete  model  found  by  reducing  the 
twentieth  order  model  (by  a  least  squares  fit  approach) ,  was 
also  unsuccessful  until  he  did  some  tuning.  In  this  study 
the  formulation  is  attempted  directly  using  the  tenth  order 
linear  state  space  model  with  the  concentrations  as  the 
states.  Since  the  control  is  carried  out  using  a  digital 
computer  the  discrete  version  of  the  control  law  is  also 
established.  In  addition  some  model  reduction  is  used  to 
simplify  the  control  action  and  to  remove  unnecessary 
measurements . 

The  emphasis  of  this  study  is  the  development  of  a 
practical  method  of  control  for  a  distillation  column.  The 
questions  to  be  answered  include  the  desirability  of  expend¬ 
ing  considerable  effort  in  developing  an  accurate  model  as 
opposed  to  an  approximate  jnodel  that  gives  some  description 

of  all  the  natural  phenomena  and  can  be  obtained  relatively 
easily,  the  importance  of  proper  control  scheme  configuration 
and  the  advantage  of  making  specific  allowance  for  on-line 
tuning . 


CHAPTER  2 


CONTROL  VARIABLE  FOR  DISTILLATE  COMPOSITION 

2.1  Introduction 

The  initial  step  in  designing  a  control  system  for  a 
process  consists  of  choosing  the  variables  that  can  be 
manipulated  to  achieve  the  desired  control.  The  system 
used  in  this  study,  a  pilot  scale  binary  distillation  column 
as  described  by  Pacey  (35)  and  Svrcek  (48)  has  five  primary 
variables  which  are  to  be  controlled.  These  are:  pressure 
(in  the  condenser),  liquid  levels  in  the  condenser  and  re¬ 
boiler  and  the  two  terminal  compositions.  The  feed  con¬ 
ditions  are  considered  to  be  the  source  of  disturbance.  The 
control  configuration  that  has  been  used  on  this  system  by 
previous  workers  (2,28,35,48)  utilizes  condenser  load  to 
control  the  pressure,  bottoms  flow  to  control  the  reboiler 
level  and  distillate  flow  to  control  the  condenser  level. 
This  leaves  the  reboiler  load  and  the  reflux  flow  to  control 
the  terminal  compositions.  This  type  of  control  system  is 
termed  energy  balance  composition  control  because  the  split 
between  overhead  and  bottoms  can  be  changed  only  indirectly 
through  a  change  in  reboiler  load  or  reflux  rate  as  the  exit 
streams  are  both  on  level  control. 

Recently  energy  balance  composition  control  has  been 
questioned.  Shinskey  (46)  states  that  since  both  product 
streams  are  on  level  control  there  is  no  direct  way  of 
manipulating  the  material  balance  and  so  better  product 
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control  could  be  obtained  by  manipulating  a  product  stream. 
Other  studies  (5,19,27,33)  have  included  tests  of  the  two 
different  systems  on  various  distillation  columns  and  shown 
that  the  control  schemes  using  a  product  flow  to  control 
composition,  termed  direct  material  balance  control  of 
composition,  gave  improved  control  behavior  for  the  systems 
tested. 

The  arguments  presented  in  favor  of  direct  material 
balance  control  are  inherent  reflux  control,  less  inter¬ 
action  and  greater  stability  for  the  controlled  system. 
Inherent  reflux  control  as  explained  by  Nisenfeld  (33) 
involves  the  reaction  of  the  system  to  a  heat  loss  distur¬ 
bance  due  to  something  like  a  sudden  decrease  in  external 
temperature.  The  system  using  the  distillate  flow  to 
regulate  the  overhead  composition  would  cause  a  decrease  in 
the  reflux  rate  (external  reflux)  but  the  total  reflux, 
internal  and  external,  would  remain  about  the  same  so  that 
the  overhead  composition  should  not  be  disturbed.  However, 
in  the  case  of  the  energy  balance  control  system,  the  reflux 
could  not  change  until  the  overhead  composition  changed  and 
thus  the  controlled  variable  would  be  disturbed. 

The  claim  of  a  decrease  in  interaction  between  the 
overhead  and  bottoms  composition  control  has  been  supported 
by  analysis  using  the  relative  gain  matrix  proposed  by 
Bristol  (5).  The  elements  in  this  matrix  are  the  ratio  of 
the  open  loop  change  in  the  controlled  variable  to  the 
change  when  the  other  controlled  variables  are  held  constant 
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under  control  for  a  change  in  the  specified  manipulated 
input.  This  was  used  to  measure  the  steady  state  inter¬ 
action  between  the  bottoms  and  overhead  composition  control 
using  alternately  the  reflux  then  the  distillate  to  control 
the  overhead  composition.  When  the  reboiler  load  is  assumed 
to  be  characterized  by  the  separation  factor  (Shinskey  (46)) 
the  relative  gain  for  the  two  alternate  overhead  composition 
control  loops  can  be  related  as  given  by  Juantoreno  and 
Romeo  ( 19 ) . 


X 


XD’R 


xd,D 


1  -  2  A 


xD’D 


2.1 


where 


X 


xd,r 


the  relative  gain  between  overhead 
composition  and  reflux  flow 


X 


xd’d 


the  relative  gain  between  overhead 
composition  and  distillate  flow 


Nisenfeld  and  Shultz  (34)  developed  an  interactive  index 


which  can  be  stated  as, 


I. 

3 


1  -  X.  . 

_ 1j_1 

X .  . 


2.2 


where 


I.  =  the  interactive  index  for  j 
3 

X .  .  =  the  relative  gain  between  i  and  j 

i,3 


Shinskey  (46)  implies  that  interaction  for  the  energy 

balance  system  is  always  worse  (except  for  X  n  =  0  or  1 

XD’ 

which  are  physically  unrealizable)  but  by  substituting 
equation  2.1  into  2.2  it  can  be  seen  that  the  interactive 
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index  is  the  same  for  both  systems.  He  also  states  that  if 

A  is  less  than  zero,  control  is  not  possible  by  the 

XD’R 

energy  balance  system.  This  is  true  but  even  for  direct 
material  balance  control  the  overhead  variable  is  best 
controlled  by  the  reboiler  load.  Thus  it  can  be  seen  that 
the  steady  state  analysis  involving  the  relative  gain  matrix 
is  not  as  useful  as  first  represented  in  choosing  the  control 
configuration  for  overhead  composition. 

One  transient  effect  which  causes  interaction  and 
degrades  control  is  the  inverse  response  undergone  by  the 
overhead  composition  for  an  increase  in  reboiler  load  which 
involves  an  initial  response  of  the  control  variable  in  the 
opposite  direction  to  the  final  steady  state  deviation. 
Rosenbrock’s  (41)  explanation  of  this  phenomenon  states  that 
the  plate  overflow  is  suddenly  increased  due  to  an  increase 
in  the  vapor  flow  which  functions  as  an  apparent  increase  in 
reflux  which  initially  causes  the  composition  in  the  column 
to  deteriorate.  The  energy  balance  control  scheme  initially 
causes  an  increase  in  the  distillate  flow  while  the  external 
reflux  remains  constant.  The  direct  material  balance  system 
immediately  increases  the  external  reflux  so  the  initial 
decrease  in  concentration  of  the  overhead  should  never  occur 
thus  eliminating  the  inverse  response. 

The  third  feature  of  the  control  system  using  the 
distillate  to  control  overhead  composition  is  the  increased 
stability  because  of  direct  control  of  the  smaller  flow. 

The  reflux  ratio  of  high  purity  product  industrial  columns 


'  w'w^  I'HI  i\  1  i  HE  1  1 


9 


is  often  over  ten.  Shinskey  (47)  states  that  in  this  case 
the  distillate  flow  is  too  small  to  control  the  level  in  the 
condenser  accumulator.  Harriot  (16)  stated  that  a  small 
error  in  the  reflux  flow  would  cause  a  much  larger  relative 
error  in  the  distillate  rate  and  thus  a  large  disturbance  to 
the  material  balance  for  the  system  with  distillate  flow  for 
level  control  however  switching  the  reflux  flow  for  level 
control  would  minimize  this  problem. 

In  view  of  the  arguments  presented  in  favor  of  direct 
material  balance  composition  control  it  is  interesting  to 
consider  why  the  energy  balance  system  has  had  wise  usage. 
This  can  be  explained  by  the  fact  that  the  reflux  is  the 
only  variable  that  can  be  used  to  directly  alter  the  vapor 
to  reflux  balance  in  the  upper  end  of  a  distillation  column 
which  in  turn  directly  changes  the  overhead  composition. 

The  alternate  system  must  rely  on  the  condenser  accumulator 
level  controller  to  change  the  reflux  flow  after  the  dis¬ 
tillate  flow  has  been  changed.  It  is  clear  that  at  least 
for  setpoint  control  the  overhead  composition  should  respond 
faster  when  the  reflux  is  directly  changed  by  the  composi¬ 
tion  control  loop. 

It  is  not  possible  to  specify  the  best  control  con¬ 
figuration  for  all  columns  because  each  different  distilla¬ 
tion  column  has  its  own  particular  characteristics.  The 
pilot  scale  distillation  column  and  the  operating  conditions 
are  such  that  they  would  diminish  the  stated  advantages  of 
the  direct  material  balance  system.  In  particular  the 


10 


reflux  ratio  is  less  than  two  so  that  the  sensitivity  of 
the  reflux  and  distillate  flows  do  not  differ  greatly.  The 
fact  that  the  column  is  located  in  the  laboratory  means  that 
no  great  change  in  heat  loss  can  be  anticipated.  Also,  the 
holdup  of  the  condenser  accumulator  is  about  the  same  as 
that  on  the  trays  so  that  the  dynamic  lag  is  small  hence 
both  reflux  and  distillate  flows  respond  at  almost  the 
same  time. 

In  view  of  the  uncertainty  of  the  virtues  of  the 
alternate  systems  the  two  systems  were  tested.  The  experi¬ 
ments  were  carried  out  with  overhead  composition  controlled 
directly  by  either  reflux  or  distillate  flow  with  the  other 
flow  maintaining  level  control  of  the  condenser  accumulator 
level.  Reboiler  load  was  maintained  constant  for  these 
tests . 

2 . 2  Experimental  Program 

The  pilot  plant  distillation  column  used  in  this  study 
has  sufficient  instrumentation  so  that  either  of  the  two 
alternate  control  schemes  for  overhead  composition  control 
as  shown  in  Figures  2.1  and  2.2  can  be  implemented  simply 
by  connecting  appropriate  pneumatic  lines. 

The  first  tests  were  carried  out  using  the  Foxboro 
capacitance  Dynalog  which  has  a  pneumatic-proportional  plus 
integral  controller  coupled  to  the  capacitance  (composition) 
measurement.  A  second  set  of  tests  were  carried  out  using 
the  DDC  (Direct  Digital  Control)  package  on  the  IBM  1800 
digital  computer  to  replace  the  controller  function  of  the 
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COOLING 


FIGURE  2.1:  Schematic  of  Overhead  Composition  Control 

Configuration  Using  Reflux  Flow  (Energy 
Balance  Control) 


COLUMN 


COOLING 


FIGURE  2.2:  Schematic  of  Overhead  Composition  Control 

Configuration  Using  Distillate  Flow  as 
Control  Variable  (Direct  Material  Balance 
Control ) 
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Dynalog.  This  was  done  to  test  the  effect  of  discrete 
control  and  to  obtain  a  more  easily  reproducible  set  of 
results.  The  reason  that  DDC  control  is  easier  to  reproduce 
stems  from  the  fact  that  input  to  the  control  function  is  in 
weight  percent  and  output  is  in  pounds  per  minute  of  flow 
rather  than  instrument  air  pressures.  Using  the  physical 
variables  gives  absolute  measurements  while  the  analog 
signals  (pressures)  are  relative  so  that  a  calibration 
change  requires  changing  the  controller  constants  to  main¬ 
tain  comparable  results. 

For  these  tests  the  distillation  column  was  subjected 
to  feed  flow  changes  and  to  a  decrease  in  feed  temperature. 
The  step  changes  in  feed  flow  were  made  to  enable  comparison 
of  performance  to  that  obtained  in  previous  control  experi¬ 
ments  (2,28,35,48).  The  feed  temperature  disturbance  was 
used  to  show  the  effect  of  an  enthalpy  change  at  the  middle 
of  the  column  in  order  to  simulate  a  change  in  heat  loss. 

Two  open  loop  responses  to  steam  flow  changes 
(changes  in  reboiler  load)  were  also  performed.  These  are 
included  to  indicate  the  degree  of  interaction  present 
between  reboiler  load  (used  for  bottom  composition  control 
in  a  latter  phase  of  this  work)  and  the  overhead  composition 
control  schemes. 

2.3  Experimental  Results 

The  initial  tests  performed  to  compare  reflux  and 
distillate  flow  as  the  manipulative  variable  with  the 
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Dynalog  controller  are  listed  in  Table  2.1.  Steady  state 
values  of  the  overhead  composition,  the  time  integral  of  the 
absolute  error  (IAE)  of  overhead  composition  along  with  the 
change  in  the  bottoms  composition  are  also  shown.  The  over¬ 
head  composition  responses  are  plotted  in  Figures  2.3  to  2.5. 
Also  included  for  comparison  are  two  runs  performed  by  Pacey. 
The  IAE  values  indicate  the  distillate  case  (distillate  flow 
used  to  control  overhead  composition)  is  generally  better. 
Figure  2.3  shows  the  effect  of  a  decrease  in  feed  flow  from 
2.46  to  1.94  pounds  per  minute.  It  can  be  seen  for  the  case 
where  reflux  flow  is  controlling  overhead  composition  that 
the  run  performed  by  Pacey  gave  poorer  control  performance 
even  though  the  two  runs  were  carried  out  with  the  same 
configuration.  In  addition  the  initial  response  of  Pacey Ts 
run  was  in  the  opposite  direction  from  the  other  two  runs 
which  would  indicate  that  the  test  started  nearer  to  the 
flooding  region.  Presumably  this  behavior  resulted  because 
the  operating  parameters  were  not  exactly  the  same  as  those 
used  by  Pacey. 

The  positive  feed  flow  step  response  depicted  in 
Figure  2.4  shows  only  small  deviation  in  the  overhead 
composition.  It  is  again  in  the  opposite  direction  from 
that  of  the  run  performed  by  Pacey. 

The  final  runs  in  which  the  analog  controller  was  used 
were  responses  using  the  alternate  overhead  composition 
control  configurations  to  a  negative  step  in  feed  temperature 
from  160°F  to  145°F.  The  response  for  the  distillate  case 


EXPERIMENTAL  RESULTS  FROM  OVERHEAD  COMPOSITION 
CONTROL  TESTS  USING  FOXBORO  DYNALOG  ANALOG  CONTROLLER 
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FIGURE  2.5:  Experimental  Response  of  Overhead  Composition 

to  a  Step  Decrease  in  Feed  Temperature  (160-145°F) 
Using  an  Analog  Controller 
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shows  some  low  amplitude  cycling.  As  the  control  valves  do 
not  have  positioners  this  was  thought  to  be  due  to  hysterisis 
in  the  reflux  flow  valve  on  level  control. 

The  tests  carried  out  using  DDC  to  perform  the  control 
function  are  listed  in  Table  2.2.  This  table  also  includes 
the  steady  state  and  IAE  of  the  overhead  composition  along 
with  the  deviation  of  the  bottoms  composition.  Table  2.3 
gives  results  from  the  run  with  the  composition  controller 
on  manual. 

The  response  of  the  overhead  composition  for  the  first 
run  using  DDC  is  shown  in  Figure  2.6.  The  disturbance  used 
was  a  feed  flow  step  from  2.4  to  2.64  pounds  per  minute. 

On  the  basis  of  maximum  deviation  in  the  open  loop  case  and 
the  IAE  in  the  closed  loop  case  it  follows  that  reflux  is 
far  superior  to  distillate  flow  for  controlling  overhead 
composition.  Comparison  of  the  controlled  response  with 
that  obtained  using  analog  control  shows  that  the  response 
is  different  in  character  which  can  again  be  attributed  to 
a  change  in  column  operating  conditions. 

The  open  loop  response  of  overhead  composition  for  a 
negative  ten  percent  feed  flow  change  shown  in  Figure  2.7 
exhibits  a  much  larger  deviation  than  occurred  for  the  step 
change  in  the  positive  direction.  This  behavior  clearly 
shows  the  nonlinearity  of  the  system.  The  fact  that  the 
closed  loop  responses  for  the  negative  step  are  unstable 
although  the  same  controller  constants  were  used  as  for  the 
positive  feed  flow  change  is  attributable  to  this  non- 
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FIGURE  2.6:  Experimental  Response  of  Overhead  Composition 

to  a  10%  Step  Increase  in  Feed  Flow  using 
DDC . 
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FIGURE  2.7:  Experimental  Response  of  Overhead  Composition 

to  a  10%  Step  Decrease  in  Feed  using  DDC. 
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linearity.  In  spite  of  the  oscillatory  response  it  is  clear 
that  the  reflux  case  has  better  control  behavior. 

Figure  2.8  shows  the  response  of  overhead  composition 
to  a  71  Fahrenheit  degree  drop  in  feed  temperature.  The 
open  loop  response  shows  that  the  overhead  composition  is 
much  more  sensitive  to  a  feed  temperature  change  when  reflux 
flow  rather  than  distillate  flow  is  the  manipulative  variable. 
This  indicates  that  distillate  flow  should  be  used  for  com¬ 
position  control  for  such  a  disturbance,  this  is  borne  out 
by  the  closed  loop  response. 

A  further  comparison  to  the  previous  work  performed  by 
Pacey  is  given  in  Figure  2.9  which  shows  responses  of  the 
overhead  composition  to  a  change  in  feed  flow  from  2.46  to 
2.92  pounds  per  minute.  It  can  be  seen  that  the  deviation 
of  the  composition  is  about  twice  as  large  when  distillate 
flow  is  the  manipulative  variable  compared  to  using  reflux. 
However,  by  tuning  the  control  the  composition  can  be  brought 
back  to  the  setpoint  in  about  the  same  period  of  time  as  for 
the  reflux  case  but  the  IAE  is  still  higher. 

Since  an  overall  control  objective  is  ideally  the 
control  of  both  overhead  and  bottoms  compositions,  it  is 
instructive  to  examine  the  deviation  in  bottoms  composition 
as  given  in  Tables  2.2  and  2.3.  The  bottoms  composition 
analysis  for  the  negative  feed  flow  steps  are  bounded  by 
zero  percent  so  they  are  not  meaningful.  It  should  be  noted 
that  the  gas  chromatograph  used  for  bottoms  composition 
determination  is  not  capable  of  detecting  concentrations 
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FIGURE  2.9:  Comparison  of  Overhead  Composition  Response 

for  an  Increase  in  Feed  Flow  from  2.46  to  2.92 
lb. /min.  for  Control  Using  Distillate  Flow 
to  Control  Using  Reflux  Flow. 
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below  about  one-half  percent  methanol  by  weight  and  results 
for  compositions  below  one  percent  are  not  reliable.  The 
deviation  for  positive  feed  flow  changes  and  feed  temp¬ 
erature  drops  indicate  that  the  direct  material  balance 
control  is  better  as  it  gives  less  interaction  to  the 
bottoms  composition.  The  deviations  in  the  bottoms  composi¬ 
tion  recorded  by  Pacey  are  inexplicably  far  less  than  com¬ 
parable  runs  shown  here. 

Also  included  in  Table  2.3  are  the  initial  and  final 
steady  states  for  the  overhead  composition  and  the  deviation 
of  the  bottoms  composition  for  open  loop  responses  to  steam 
flow  changes  (reboiler  load  change)  when  distillate  and  re¬ 
flux  flow  are  manipulated  variables.  The  results  indicate 
that  the  overhead  is  least  affected  by  changes  in  steam  rate 
when  distillate  flow  is  used  for  composition  control  rather 
than  reflux  flow. 

2.4  Discussion  of  the  Experimental  Results 

The  results  of  the  single  variable  control  scheme  study 
show  evidence  that  if  only  the  overhead  composition  is  to  be 
controlled  it  would  be  better  to  use  the  reflux  flow.  The 
runs  using  the  analog  controller  generally  indicate  some 
improvement  for  the  case  using  distillate  flow  to  control 
the  overhead  composition  but  the  improvement  is  not  great 
enough  to  be  conclusive.  If,  however,  the  DDC  runs  for  feed 
flow  steps  are  examined  it  can  be  seen  that  the  difference 
is  quite  pronounced  in  favor  of  the  reflux  case.  Since  feed 
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flow  steps  are  the  basic  disturbance  used  on  the  column  it 
would  indicate  that  the  reflux  flow  is  the  variable  that 
should  be  used  in  order  to  get  the  best  control. 

The  expansion  of  the  problem  to  include  control  of 
both  terminal  compositions  changes  this  conclusion.  The 
fact  that  there  is  less  interaction  in  both  directions  when 
distillate  flow  is  used  is  strongly  in  favor  of  this  con¬ 
figuration.  It  is  because  of  this  fact  that  the  direct 
material  balance  control  scheme  utilizing  distillate  flow 
and  reboiler  load  as  manipulated  variables  to  control  over¬ 
head  and  bottoms  compositions  respectively  was  selected  for 
use  in  the  subsequent  control  study  concerned  with  terminal 
composition  control. 


. 


CHAPTER  3 


MODELLING  THE  DISTILLATION  COLUMN 

3.1  Introduction 

Some  attempts  have  been  made  to  model  the  University 
of  Alberta  pilot  plant  distillation  column  with  varying 
success.  Svrcek  (48)  derived  a  relatively  complex  model 
from  theoretical  considerations  and  experimental  data.  He 
was  able  to  predict  the  experimental  response  of  the  column 
with  reasonable  accuracy.  The  studies  that  followed  were 
concerned  with  fitting  transfer  function  models  to  pulse  test 
data.  Berry  (2)  and  Pacey  (35)  used  first  order  plus  time 
delay  transfer  function  forms.  They  found  that  the  system 
was  non-linear  to  the  extent  that  the  models  were  only 
accurate  for  small  excursions  from  the  steady  state. 

Chanh  (7)  fitted  second  order  transfer  functions  with  a  time 
delay  to  pulse  test  data  and  again  found  the  system  quite 
non-linear . 

McGinnis  (28)  modelled  the  column  using  the  linearized 
liquid  enthalpy  relationship  and  a  first  order  lag  to  re¬ 
present  the  liquid  dynamics  for  each  tray.  His  assumption 
that  the  vapor  flow  change  was  negligible  is  questionable  as 
he  used  reboiler  load  changes  as  one  of  his  control  variables 
in  the  subsequent  control  study.  He  overcame  the  problem  by 
adding  experimental  values  of  the  condenser  load  change  to 
his  runs  and  so  obtained  fairly  good  agreement  between  the 
predicted  response  from  the  model  and  the  experimental 
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response  of  the  equipment.  Another  shortcoming  of  this 
model  was  that  it  gives  good  results  only  for  conditions 
associated  with  high  purity  bottoms  composition. 

In  view  of  the  shortcomings  of  previous  studies  this 
work  on  modelling  was  undertaken  to  examine  the  system  and 
the  assumptions  and  to  attempt  to  find  a  simple  model  with 
a  wide  range  of  use.  The  first  item  examined  was  whether 
liquid  enthalpies  or  liquid  compositions  should  be  used  to 
characterize  the  stages  (states).  Levy  et  al  (22)  applied 
modal  analysis  to  the  models  resulting  from  using  the 
energy  and  material  balances  for  a  binary  distillation 
column  separating  an  acetone-benzene  mixture.  They  showed 
that  the  component  mass  balances  were  the  best  to  use  to 
characterize  the  system.  Although  this  analysis  is  somewhat 
dependent  on  the  physical  system  it  is  applicable  to  the 
University  of  Alberta  distillation  column. 

In  order  to  evaluate  different  models  it  is  necessary 
to  calculate  the  transient  response  of  the  states  to  some 
disturbance  input.  McClune  and  Gallier  (27)  who  used  the 
IBM  program  CSMP/360  (Continuous  System  Modelling  Program/ 360) 
in  designing  a  distillation  column  control  system  have 
discussed  the  advantages  of  using  a  digital  simulation 
program.  An  improved  version  of  this  program  CSMP  III  (57) 
was  used  in  this  study  to  calculate  the  predicted  responses 
from  the  models  for  comparison  to  the  experimental  data. 

Modelling  of  distillation  processes  has  received  much 
attention  in  the  past  such  as  the  articles  by  Rosenbrock 
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(42,43,44).  The  literature  dealing  with  modelling  with 
respect  to  distillation  has  been  categorized  by  McGinnis  and 
by  Levy  et  al  showing  the  types  of  models  used  and  the 
assumptions  made.  The  reports  are  generally  case  studies 
because  distillation  processes  cannot  be  generallized  to  any 
great  extent.  For  that  reason  the  emphasis  in  this  study 
will  be  on  obtaining  a  model  for  the  University  of  Alberta 
distillation  column,  then  the  model  can  be  examined  for  its 
general  applicability. 

3.2  Description  of  the  Experimental  Equipment 

The  University  of  Alberta  distillation  column  is  nine 
inches  in  diameter  with  glass  walls  to  enable  viewing  of  the 
process.  It  has  eight  bubble  cap  trays,  a  thermosyphon 
reboiler  and  a  total  condenser.  A  detailed  description  of 
the  distillation  column  and  the  associated  equipment  is  given 
by  Svrcek  (48)  and  by  Pacey  (35). 

The  column  is  operated  with  a  feedrate  of  about  2.4 
pounds  per  minute  of  methanol-water  mixture  containing  47 
weight  percent  methanol.  This  is  separated  to  about  2  per¬ 
cent  in  the  bottoms  and  97  percent  in  the  overhead.  A 
complete  set  of  typical  operating  conditions  are  given  in 
Appendix  A.  The  overall  control  objective  is  to  maintain 
both  the  overhead  and  bottoms  compositions  at  some  constant 
value . 

The  column  has  extensive  instrumentation  to  enable  a 
wide  range  of  control  studies  to  be  undertaken.  Temperature 
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measurements  are  taken  on  all  the  trays  and  in  the  condenser 
and  reboiler  in  the  liquid  phases.  The  feed,  reflux,  bottoms 
and  distillate  flowrates  are  measured  as  are  their  temp¬ 
eratures.  Inlet  and  outlet  temperatures  and  the  flow  of  the 
steam  and  condenser  cooling  water  are  measured.  All  these 
measurements  as  well  as  the  liquid  compositions  in  the  re¬ 
boiler  and  condenser  are  transmitted  to  the  IBM  1800  digital 
computer.  The  bottoms  composition  is  measured  by  an  indus¬ 
trial  gas  chromatograph  controlled  by  the  IBM  1800  on  a  cycle 
of  approximately  three  minutes . 

The  column  is  outfitted  with  a  complete  set  of  analog 
controllers  to  control  the  flowrates  as  well  as  reboiler  and 
condenser  levels  and  the  feed  and  reflux  temperatures.  The 
setpoints  of  the  principle  pneumatic  flow  controllers  can  be 
set  from  the  IBM  1800  in  a  supervisory  control  mode. 

A  schematic  diagram  of  the  column  giving  details  of 
the  controls  and  measurements  is  given  in  Appendix  A. 

When  the  outputs  from  the  computer  are  maintained 
constant  the  column  control  is  said  to  be  open  loop  even 
though  the  local  variables  such  as  level  are  controlled 
under  analog  controllers.  A  simple  schematic  of  the  column 
configuration  is  shown  in  Figure  3.1 

3.3  Derivation  of  the  Models 

The  usual  approach  to  modelling  a  distillation  column 
begins  by  considering  the  column  to  consist  of  a  number  of 
stages  connected  by  mass  and  energy  flows.  Figure  3.2, 
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FIGURE  3.1:  Schematic  of  the  Distillation  Column  Control  Scheme 
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Representation  of  the 


Stage . 


FIGURE  3.2: 


32 


which  is  a  representation  of  a  typical  stage,  shows  the 
variables  that  will  be  used.  The  stage,  used  to  represent 
a  tray,  reboiler  or  condenser  will  be  considered  to  have  two 
distinct  phases,  liquid  and  vapor.  Each  phase  is  assumed  to 
be  perfectly  mixed.  In  addition  other  assumptions  can  be 
made  about  the  stage  such  as  a  fixed  relation  between  the 
physical  properties  of  the  liquid  phase  and  those  of  the 
vapor  phase. 

For  each  stage  material  and  energy  balances  can  be 
written.  For  the  time  varying  case  the  material  balances 
can  be  stated  as : 


3.1 


3.2 


The  energy  balance  is : 


3.3 


Since  the  column  contains  eight  trays  plus  reboiler 
and  condenser  there  are  ten  stages  in  all.  These  balances 
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yield  thirty  equations  which  contain  140  unknowns  (fourteen 
per  stage).  In  order  to  reduce  the  number  of  unknowns  to 
thirty  so  that  the  equations  can  be  solved  some  assumptions 
must  be  made. 

Certain  variables  are  classified  as  known  inputs 
because  they  are  measured  or  easily  calculated  and  directly 
externally  controllable.  Included  in  this  group  are  feed 
flow,  enthalpy  and  composition,  reboiler  load,  condenser 
load,  reflux  flow,  distillate  flow  and  bottoms  flow.  Thus 
the  feed  and  sidestream  drawoffs  are  known  for  each  stage 
reducing  the  number  of  unknowns  by  fifty  to  ninety.  If  some 
assumption  is  made  about  heat  loss  for  the  stages  such  as 
constant  heat  loss  the  heat  transfer  for  each  stage  is  known 
leaving  eighty  unknowns. 

It  was  stated  previously  that  physical  properties  of 
the  vapor  could  be  assumed  to  have  some  functional  relation¬ 
ship  with  those  of  the  liquid  phase.  This  means  that  vapor 
composition  and  enthalpy  could  be  found  from  the  liquid 
composition  and  enthalpy  respectively  thus  reducing  the 
unknowns  by  twenty  to  sixty. 

Other  assumptions  can  be  made  with  respect  to  liquid 
and  vapor  dynamics.  For  example,  assuming  the  vapor  dynamics 
negligible  reduces  the  unknowns  to  fifty  because  vapor  hold¬ 
ups  can  be  assumed  to  be  the  steady  state  values.  On  the 
other  hand,  liquid  dynamics  can  be  approximated  by  using 
weir  equations  or  by  assuming  constant  liquid 

This  gives  a  relationship  between  holdup  and 


holdup . 
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liquid  flow  and  reduces  the  unknowns  by  ten  to  forty. 

The  last  ten  unknowns  may  be  eliminated  by  assuming 
some  relationship  for  the  vapor  flowrate  based  on  known 
inputs  such  as  the  reboiler  load. 

Even  though  the  equations  can  now  be  solved  (providing 
all  the  functional  relationships  assumed  can  be  found)  there 
is  still  the  problem  that  on  the  actual  column  some  variables 
such  as  liquid  holdup  are  not  measured.  Control  schemes  such 
as  optimal  control  require  that  all  the  " states"  be  available 
(measured).  States  in  this  case  are  the  unknowns  remaining 
in  the  stage  equations.  It  might  be  possible  to  estimate 
some  states  by  an  estimation  technique  but  in  the  case 
of  liquid  holdups  it  is  doubtful  if  such  a  procedure 
would  prove  satisfactory.  One  way  to  eliminate  the 

unwanted  variables  is  to  make  more  empirical  assumptions. 
Another  approach  is  to  use  numerical  model  reduction 
techniques  such  as  Marshall’s  (26)  modal  reduction  method  or 
Anderson’s  (1)  least  squares  reduction  technique.  Both 
approaches  to  eliminating  states  have  disadvantages.  The 
empirical  approach  is  based  on  intuition  and  there 
criteria  to  determine  whether  one  assumption  is  better  than 
another  and  no  guide  as  to  what  assumptions  can  justifiably 
be  made.  The  modal  and  least  squares  reduction  techniques 
require  a  linear  system  which  generally  means  the  system 
equations  must  be  linearized.  In  addition  there  is  the 
problem  of  choosing  which  states  to  retain. 
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Since  the  column  is  operated  under  pressure  control 
via  the  condenser  load  it  can  be  assumed  that  the  pressure 
at  each  stage  remains  constant.  As  a  result  of  this 
assumption  it  follows  that  if  temperature  (or  liquid 
enthalpy)  of  a  stage  is  known  then  the  composition  of  the 
liquid  phase  can  be  established  or  knowing  the  composition 
the  enthalpy  can  be  calculated.  Thus  retaining  both 
enthalpies  and  compositions  as  states  is  meaningless  because 
they  are  related.  The  problem  is  in  choosing  which  variable 
to  retain.  As  no  clear  choice  exists  each  was  used  to  derive 
a  model  of  the  distillation  column.  In  the  derivations  the 
liquid  flow  dynamics  were  neglected  in  order  to  simplify  the 
model  so  that  reduction  was  not  necessary.  Equations  for 
liquid  dynamics  are  presented  for  completeness. 

3.3.1  Enthalpy  Model 

The  simplest  dynamic  model  that  can  predict  the  response 
of  all  ten  stages  must  be  of  tenth  order.  Considering  the 
pressure  as  constant  and  the  liquid  phase  of  each  stage  to 
be  at  the  bubble  point  it  is  possible  to  specify  the  column 
compositions  from  temperatures  or  enthalpies.  Possible 
shortcomings  in  this  scheme  include  unaccounted  for  pressure 
deviations,  non-saturated  conditions,  variable  heat  loss  and 
incomplete  or  unreliable  enthalpy  data.  Direct  manipulation 
of  the  heat  balance  to  control  the  primary  variables  namely 
the  terminal  compositions  is  possible  using  this  model. 

A  summary  of  the  various  assumptions  necessary  to  form 
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a  model  containing,  as  the  only  unknowns,  the  liquid  enthal¬ 
pies  on  each  stage  are: 


1.  Constant  pressure  =  constant 

implies  x.  =  f, (h.)  and  h.  =  f0(T.) 
r  ill  1  2  1 

2.  Negligible  liquid  dynamics  =  0 

3.  Negligible  vapor  dynamics  =  0 

4.  Constant  heat  loss  q^  =  0 

5.  Physical  properties  of  the  vapor  are 
related  to  those  of  the  liquid  Y.  =  f~(x.) 

1  o  1 


where 


f^,  ^2’  ^3  an<^  ^4  <^eno"t:e  functional 
relationships  and  i  =  1,2, . ,10. 


In  addition  to  the  above  the  vapor  flow  can  be  esti¬ 
mated  from  the  heat  input  on  the  assumption  that  heat  is 
transferred  instantly  up  the  column  by  changes  in  vapor  flow. 
McGinnis  (28)  used  the  assumption  that  the  vapor  flowrate 
remained  constant  which  is  reasonable  so  long  as  the  re¬ 
boiler  load  is  kept  constant.  However,  as  the  reboiler  load 
is  one  of  the  primary  control  variables  it  is  considered  a 
poor  assumption.  The  actual  relationship  used  for  calcula¬ 
ting  the  vapor  rate  from  any  stage  includes  the  effect  of 
heat  losses  from  each  tray  and  the  effect  of  a  non-saturated 
feed  giving: 


[Qr  +  ?  Oj  +  f. 

3=1 


Ah 


F. 

3 


hP 


3.4 
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where  QR  -  reboiler  load  =  Qq  “  Gq 

-  heat  loss  from  each  tray 
Ah-p  -  the  difference  between  the  feed 

Jl  • 

3 

enthalpy  and  the  saturated  feed  enthalpy. 

In  order  to  simplify  the  resulting  equations  a  relation 
between  the  liquid  and  vapor  enthalpy  is  needed.  McGinnis 
used  a  linear  relation  based  on  the  equilibrium  data  for 
the  methanol-water  system.  His  model  was  in  linear  per¬ 
turbation  form  justifying  his  use  of  a  linear  relation. 

The  intent  in  this  work  is  to  develop  a  non-linear  model 
which  has  a  wider  range  of  valid  operating  conditions.  The 
method  chosen  consisted  of  using  pseudo-equilibrium  data 
taken  from  steady  state  runs  performed  by  McGinnis  (28)  and 
Svrcek  (48). 


The  time  varying  enthalpy  balance  for  a  general  stage 


is  given  by: 


h .  WT  +  h.  W,  +  H.  W„ .  +  H.  W 

i  L.  l  L.  i  Vi  i  Vi 
i  l 

=  Vi-1  Hi-1  +  Li+1  hi+l  +  Fi  hFi 


3.5 


-  (Li  +  SLi)  hi  (Vi  +  Syi)  Hi  +  Qi 


The  terms  containing  the  derivatives  of  the  holdups 
can  be  neglected  by  considering  liquid  and 

vapor  holdup  to  be  constant.  The  term  containing  the  vapor 
holdup  was  retained  in  spite  of  the  fact  that  the  magnitude 
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of  the  vapor  holdup  is  small  because  the  derivative  of  the 
vapor  enthalpy  is  relatively  large.  In  order  to  retain  that 
term  the  relationship  between  the  derivative  of  the  liquid 
and  vapor  enthalpy  was  linearized 


where 


p=constant 


3.6 


This  can  be  justified  by  the  fact  that  the  small  vapor 
holdup  reduces  the  relative  non-linear  effect  of  this  term 
and  by  the  fact  that  this  effects  only  the  time  response  and 
not  the  steady  state  values. 

The  enthalpy  balance  can  now  be  written  as: 


[V .  ,  H .  +  L .  ,  h .  ,,  +  F.h-,-,. 

l-l  l-l  i+l  l+l  l  Fi 


-  (L.  +  ST.)  h.  -  (V.  +  S„.)  H.  +  Q.] 
l  Li  l  l  Vi  l  xi 


CWLi  +  Bli  WVi] 
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This  equation  for  a  general  stage  can  now  be  modified 
to  fit  each  of  the  stages  in  the  actual  column.  This  is 
carried  out  in  Appendix  B. 


3.3.2  Concentration  Model 

The  concentration  model  utilizes  the  liquid  composition 
(concentration  of  methanol)  to  characterize  stage  behavior. 
This  model  has  the  advantage  of  explicitly  using  the  terminal 
compositions  as  states. 

The  assumptions  used  for  this  model  can  be  summarized 
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as  follows: 


1. 


2. 

3. 

4. 


5. 


Constant  pressure  =  constant 

implying  xi  =  fgClP),  Hi  =  f^CY^ 

Negligible  liquid  dynamics  = 

Negligible  vapor  dynamics  = 

Constant  heat  loss  q.  =  constant 
Vapor  composition  is  a  function  of 


liquid  composition  Y.  =  f„(x.) 

1  o  1 


o 

0 


The  same  expression  for  vapor  flow  is  used  for  the 
concentration  model  as  was  used  for  the  enthalpy  model. 

The  functional  relationship  between  vapor  and  liquid 
compositions  f ^ ,  was  found  in  the  same  manner  as  the 
relation  between  liquid  and  vapor  enthalpy.  The  steady 
state  data  from  Svrcek  (48)  was  smoothed  to  form  a  pseudo¬ 
equilibrium  curve  and  represented  by  sufficient  points 
to  enable  linear  interpolation  for  simulation  purposes. 

The  component  balance  can  be  written  as: 


WT  .  x.  +  Ww .  Y.  +  WT  .  x.  +  W„ •  Y. 
Li  l  Vi  l  Li  l  Vi  l 


=  L.,-.  x.,n  +  V.  ,  Y.  ,  -  (L.  +  ST  .  )  x. 
l+l  l+l  l-l  l-l  l  Li  l 


3.8 


-  (V.  +  S„. )  Y.  +  F.  xr. 
l  Vi  l  l  Fi 


The  vapor  holdups  were  considered  negligible  as  they 
were  less  than  a  tenth  of  the  liquid  holdups  so  with  some 
rearrangment  equation  3.8  becomes: 
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Li 


3.9 


The  modification  of  this  equation  to  fit  the  actual  stages 
is  carried  out  in  Appendix  B. 

3.4  Linearization  of  the  Concentration  Model 

The  control  law  chosen  for  this  study  is  restricted  to 
linear  systems  so  the  model  was  linearized  by  taking  the 
linear  terms  of  the  Taylor  Series  expansion.  This  is 
necessary  because  the  component  balances  contain  product 
terms . 

In  addition  it  is  necessary  to  have  linear  relation¬ 
ships  for  vapor  flows  and  compositions.  These  relations 
are  given  in  Appendix  B. 

The  resulting  equations  for  the  system  which  appear  in 
Appendix  B  are  of  the  form 


f (z ,u ,d) 


3.10 


where 


u 


z 


a  vector  of  the  states  [x^ ,X£ , • . • ,x^q] 
a  vector  of  control  inputs  [D,Q^] 


5  •  •  •  5 


d 


a  vector  of  disturbance  inputs  [F,Xp,hp] 


This  can  be  put  in  the  form: 


z 


Az  +  Bu  +  Cd 


3.11 


where 


A,  B  and  C  -  coefficient  matrices 
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The  coefficient  matrices  are  constants  for  a  particular 
initial  steady  state.  A  set  of  values  for  a  typical  steady 
state  are  given  in  Appendix  A. 

3 . 5  Results 

The  test  program  involved  a  series  of  experimental 
runs  for  different  disturbances  with  the  results  being 
compared  to  the  predicted  response  from  the  models.  The 
models  were  tested  for  disturbances  in  feed  flow,  feed 
temperature,  reboiler  load  and  distillate  flow. 

The  linear  enthalpy  model  was  tested  first.  The 
results  from  this  were  very  poor.  This  was  attributed  to 
the  high  sensitivity  of  the  enthalpy  of  the  saturated 
solution  versus  concentration  relationship.  A  plot  of  this 
relationship  is  given  in  Appendix  A.  This  means  that  a  very 
small  error  in  the  heat  balance  will  make  a  large  difference 
in  the  concentration  at  the  higher  methanol  concentrations. 

No  results  from  this  model  have  been  presented  as  the  results 
were  quite  poor.  This  model  was  abandoned  in  favor  of  the 
concentration  model. 

The  non-linear  concentration  model  was  next  tested  for 
its  ability  to  accurately  predict  the  response  of  the  distil¬ 
lation  column.  The  initial  and  final  steady  state  values  can 
be  compared  by  examining  the  experimental  values  given  in 
Table  3.1  and  the  values  calculated  from  the  model  given  in 
Table  3.2.  Although  these  results  were  far  better  than  those 
from  the  enthalpy  model  they  are  not  as  good  as  the  various 
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TABLE  3,2 

PREDICTED  FINAL  STEADY  STATE  VALUES  FROM  TENTH  ORDER  NONLINEAR 
CONCENTRATION  MODEL  FOR  OPEN  LOOP  RESPONSES 
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results  from  the  models  presented  by  previous  workers.  This 
model  is  not  as  dependent  on  the  accuracy  of  the  input  con¬ 
ditions  as  was  the  case  of  those  of  previous  workers.  The 
previous  models  were  based  on  a  specific  set  of  operating 
conditions . 

The  linearized  versions  of  this  model  were  taken 
from  the  steady  states  calculated  from  the  non-linear  model 
using  the  input  conditions  of  the  experimental  runs.  The 
first  three  linear  models  were  based  on  the  initial  con¬ 
ditions  of  the  ten  percent  step  increase  in  reboiler  load, 
the  five  percent  increase  in  distillate  flow  and  the  ten 
percent  increase  in  feed  flow.  The  final  linearization  was 
done  around  the  final  steady  state  of  the  ten  percent  step 
increase  in  feed  flow.  The  steady  state  conditions  used 
for  linearization  of  the  models  are  shown  in  Table  3.3.  From 
Figure  3.3,  which  shows  the  predicted  responses  from  the 
linearized  models  and  the  experimental  data  for  a  ten  percent 
step  disturbance  in  feed  flow,  it  can  be  seen  that  the  pre¬ 
dicted  response  is  highly  dependent  on  the  steady  state 
around  which  the  linearization  is  performed.  Although  this 
indicates  that  the  system  cannot  be  represented  by  a  linear 
model  for  any  sizeable  deviation  from  the  original  steady 
state  the  linear  model  was  used  in  the  subsequent  control 
study.  The  reason  for  this  is  simply  that  no  method  of  in¬ 
corporating  the  nonlinearities  of  the  system  in  the  type  of 
control  scheme  being  considered  was  found.  If  the  system 
can  be  kept  under  tight  control  the  deviation  from  the 
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TABLE  3 . 3 

STEADY  STATE  CONDITIONS  USED  FOR  LINEARIZATION 


MODEL 

1 

2 

3 

4 

Feed  Flow 
( lb/min) 

2.4 

2.4 

2 . 4 

2.64 

Reboiler  Duty 
( BTU/min) 

1782.1 

1782.1 

1783.0 

1783.0 

Distillate  Flow 
( lb/min) 

1. 092 

1.  058 

1.1 

1.1 

Feed  Composition 
(Wt.  Pet.  MEOH) 

46.74 

46 . 31 

44.86 

44.86 

Feed  Temperature 
C°  F) 

161.0 

161.  0 

161.0 

161.0 
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initial  steady  state  would  be  small  so  the  linear  model 
should  be  adequate. 

The  disturbance  used  to  produce  the  response  shown 
in  Figure  3.3  was  a  ten  percent  increase  in  feed  flow.  The 
non-linear  model  gives  the  best  prediction  of  the  bottoms 
response.  The  overhead  response  is  comparatively  small  at 
these  conditions  because  the  overhead  composition  is  near 
the  pinch  region.  The  simulation  predicts  an  even  smaller 
overhead  response.  Examining  the  response  reveals  that  the 
process  has  a  time  delay  of  about  ten  minutes.  The  actual 
system  time  constants  lie  between  those  of  the  linearized 
models  but  seem  to  be  longer  than  that  of  the  non-linear 
model . 

The  response  for  a  negative  change  in  feed  flow  is 
not  plotted  as  the  composition  of  both  the  overhead  and 
bottoms  streams  went  out  of  the  range  of  the  measuring 
systems . 

The  response  of  the  system  to  a  ten  percent  increase 
in  steam  flow  (change  in  reboiler  load)  as  shown  in  Figure 
3.4  does  not  agree  well  with  the  predicted  response  from  the 
concentration  model.  The  overhead  composition  changed  in 
the  opposite  direction  from  that  predicted.  From  the 
response  to  a  negative  ten  percent  change  in  steam  as  shown 
in  Figure  3.5  it  can  be  seen  that  the  change  in  overhead 
composition  is  in  the  same  direction  as  for  the  positive 
step  disturbance.  No  linear  model  can  predict  such  behavior 
as  the  principle  of  superposition  must  hold.  The  physical 
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0  40  80  120 
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Comparison  of  Simulated  Responses  Predicted  by 
Concentration  Models  and  Experimental  Response 
for  a  Ten  Percent  Step  in  Feed  Flow 
(2.4-2.64  LB/MIN) 


FIGURE  3.3: 
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OVERHEAD 


TIME  IN  MINUTES 

FIGURE  3.4:  Comparison  of  Simulated  Response  Calculated 

From  Linear  Concentration  Model  to  Experimental 
Response  for  a  Ten  Percent  Step  in  Steam  Flow 
to  the  Reboiler  (1.85-2.035  LB/MIN) 


COMPOSITION  CHANGE  (wtpct) 


49 


TIME  IN  MINUTES 


Comparison  of  Response  Predicted  by  Linear 
Concentration  Model  to  Experimental  Response 
for  a  Ten  Percent  Step  Decrease  in  Steam  Flow 
(1.85-1.665  LB/MIN) 


FIGURE  3.5: 
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reason  for  this  was  deemed  to  be  the  change  in  efficiency 
due  to  the  holdup  change  that  accompanies  a  change  in  steam 
flow.  The  negative  change  in  reboiler  load  gives  better 
agreement .  In  comparing  this  response  to  that  from  the  feed 
disturbance  it  can  be  seen  that  the  time  delay  is  no  longer 
present  and  the  response  is  fairly  close  in  time  for  the 
simulation  and  the  experimental  results.  However ,  the  mag¬ 
nitude  of  the  responses  differ  by  a  factor  of  two. 

The  step  change  disturbance  of  ten  percent  in  the 
distillate  flow  resulted  in  the  composition  measurement 
instruments  going  off  scale.  A  five  percent  increase  was 
then  carried  out.  The  response  of  the  terminal  compositions 
along  with  the  response  predicted  by  the  non-linear  and 
linear  models  are  shown  in  Figure  3.6.  The  final  steady 
state  agreement  is  fair  but  like  the  response  to  a  feed  flow 
step  there  is  a  time  delay  that  has  not  been  accounted  for 
in  the  model.  It  is  possible  that  this  delay  is  due  mainly 
to  the  liquid  dynamic  lags  in  the  liquid  flows  on  the  trays. 
The  negative  change  of  ten  percent  in  the  distillate  flow 
resulted  in  the  response  shown  in  Figure  3.7.  The  effect  of 
an  increase  in  plate  holdup  up  to  the  point  where  flooding 
starts  is  evident  here.  The  models  predict  concentrations 
that  are  lower  for  both  overhead  and  bottoms  up  to  flooding 
then  overhead  composition  drops  off. 

The  final  disturbance  is  a  change  in  feed  temperature 
from  161° F  to  142° F  as  shown  in  Figure  3.8.  The  response  of 
the  system  is  very  small  but  appears  to  be  in  the  correct 
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TIME  IN  MINUTES 


FIGURE  3.6:  Comparison  of  Response  Predicted  by  Linear 

Concentration  Itodel  to  Experimental  Response 
for  a  Five  Percent  Step  in  Distillate 
(1.058-1.108) 
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o  40  80  120 


TIME  IN  MINUTES 

FIGURE  3.7:  Comparison  of  Simulated  Responses  Calculated 
From  Concentration  Models  to  Experimental 
Response  for  a  Ten  Percent  Step  Decrease  in 
Distillate  Flow  (1.1-.99  LB/MIN) 


COMPOSITION  CHANGE  (  wt  pcf ) 
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TIME  IN  MINUTES 


FIGURE  3.8:  Comparison  of  Response  Predicted  by  the 

Nonlinear  Concentration  Model  to  the  Experimental 
Response  for  a  Step  Decrease  in  Feed  Temperature 
(161-142°F) 
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direction.  This  shows  the  relative  insensitivity  of  the 
system  to  mid-column  energy  changes. 

It  should  be  noted  that  the  high  noise  level  in  the 

bottoms  composition  measurement  necessitated  some  smoothing 

of  the  plots.  This  was  done  ”by  eye  "  on  the  unfiltered 

measurements.  Improvement  to  the  gas  chromatograph  measuring 

system  or  heavy  filtering  is  necessary  in  order  to  obtain 

reasonable  control.  The  predicted  deviations  are  given  in 

Table  3.4. 

3 . 6  Discussion 

The  predictions  of  the  model  used  here  are  not  as 
satisfactory  as  those  obtained  previously  by  Svrcek,  Pacey , 
Berry,  Chanh  and  McGinnis.  This  means  that  agreement  could 
be  greatly  improved  by  using  more  experimental  data  and  by 
adding  greater  sophistication  to  the  model  by  including  such 
effects  as  the  liquid  dynamic  lags  and  perhaps  use  an 
efficiency  as  a  function  of  loading  rather  than  the  pseudo¬ 
equilibrium  curve.  The  question  still  remains  as  to  how 
much  effort  needs  to  be  expended  in  model  development  in 
order  to  enable  reasonable  design  of  the  control  system. 

The  answer  to  this  question  depends  strongly  on  the  design 
technique  used  to  find  a  control  scheme.  All  the  previous 
workers  found  it  necessary  to  tune  the  resulting  control 
parameters  when  the  scheme  was  implemented  on  the  experi¬ 
mental  equipment. 

The  fact  that  fairly  large  disturbances  were  used  on  a 
highly  non-linear  piece  of  equipment  implies  that  the  linear 
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model  cannot  be  expected  to  give  good  results.  In  view  of 
the  phenomena  that  were  assumed  negligible  the  model  gives 
reasonable  representation  of  system  behavior. 

From  a  practical  point  of  view  the  physical  features 
such  as  measurement  noise,  accuracy  and  control  action 
reliability  could  be  more  important  than  model  accuracy  in 
the  final  design  steps.  The  measurement  noise  associated 
with  the  gas  chromatograph  which  measures  the  bottoms 
composition  is  a  significant  shortcoming  of  the  experimental 
system.  This  is  one  of  the  primary  variables  and  the 
recorded  measurement  exhibits  a  noise  level  up  to  100  percent 
or  more  of  the  actual  measurement  at  the  low  compositions 
being  measured.  Other  measurement  problems  include  the 
effect  of  density  changes  on  flow  readings  measured  by  a 
pressure  drop  across  quadrant  and  square-edged  orifices  and 
the  noise  from  the  temperature  measurements.  In  addition, 
control  accuracy  is  affected  by  the  fact  that  the  control 
valves  do  not  have  positioners. 

It  appears  that  the  assumptions  made  were  at  least 
partially  effective.  The  vapor  flowrate  is  somewhat  affected 
by  heat  loss  changes  in  the  column.  Since  the  insulating 
plexiglass  was  removed  some  variation  here  could  be  expected. 
The  idea  of  pseudo-equilibrium  has  the  advantage  of  enabling 
a  wide  range  of  operating  conditions  but  it  does  not  account 
for  the  change  in  efficiencies  throughout  the  column  affected 
by  different  loading  so  is  limited  in  accuracy.  The  fact 
that  the  conditions  used  for  the  linearization  were  taken 


< 


57 


from  the  model  rather  than  from  experimental  values  leads  to 
error  but  does  give  considerable  ease  in  changing  conditions. 
It  would  be  possible  to  continuously  calculate  the  controller 
based  on  the  updated  model  using  this  method  but  for  this 
approach  to  be  advantageous  the  accuracy  of  the  non-linear 
model  would  have  to  be  improved. 


CHAPTER  4 


EXPERIMENTAL  EVALUATION  OF  A  CONTROL  SCHEME 

4.1  Introduction 

For  a  distillation  column  it  is  possible  to  control 
the  concentration  of  a  component  in  both  the  overhead  and 
bottoms  streams.  However,  in  practice,  control  is  made 
difficult  by  the  interaction  present  in  the  system  and  by 
sluggish  response  of  the  terminal  compositions.  The  intent 
here  is  to  present  a  possible  solution  to  these  problems  by 
a  modification  of  the  control  method  proposed  by  Greenfield 
and  Ward  (13,14)  and  to  apply  the  method  using  the  conc¬ 
entration  model  derived  in  the  previous  chapter.  The  control 
behavior  was  first  studied  by  simulation  using  the  model 
before  testing  it  on  the  actual  equipment  for  experimental 
verification . 

Previous  workers  at  the  University  of  Alberta  have 
tested  several  schemes  on  the  distillation  column.  Svrcek 
(48),  who  was  concerned  mainly  with  the  dynamics  of  the 
system,  tried  only  overhead  composition  control  with  a  three 
mode  feedback  controller.  Pacey  (35)  also  controlled  only 
the  overhead  composition  but  he  showed  that  considerable 
improvement  could  be  made  by  adding  feedforward  action.  He 
made  use  of  a  method  described  by  Luyben  (25)  to  find  feed¬ 
forward  control  functions.  Chanh  (7)  studied  overhead  com¬ 
position  control  using  a  measurement  of  an  intermediate  tray 
temperature  between  the  feed  tray  and  the  condenser.  Use  of 
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an  intermediate  tray  temperature,  however,  leads  to  final 
steady  state  offset. 

The  simultaneous  control  of  both  bottoms  and  overhead 
compositions  was  implemented  by  Berry  (2).  He  showed  that 
for  two  point  composition  control  the  interaction  between 
the  terminal  compositions  and  the  manipulated  variables, 
reboiler  load  and  reflux  flow,  degraded  the  control  per¬ 
formance.  He  then  tested  two  schemes  to  remove  or  diminish 
the  interaction.  The  first  proposed  by  Rijnsdorp  and  Van 
Kampen  (40),  used  a  ratio  controller  to  control  the  reflux 
ratio  while  the  ratio  setpoint  was  adjusted  by  the  overhead 
composition  controller.  The  second  scheme,  non-interacting 
control,  involved  the  addition  of  decoupling  elements  as 
outlined  by  Zalkind  (55,56)  and  by  Luyben  (24).  Unlike  the 
first  scheme,  the  second  is  based  on  the  use  of  a  transfer 
function  model.  The  results  indicate  that  the  simpler 
method  (reflux  ratio  scheme)  gave  results  comparable  to  the 
non-interacting  scheme  but  both  schemes  gave  improvement 
over  the  case  where  there  was  no  compensation  for  the 
interaction. 

McGinnis  (28), using  a  theoretical  model,  implemented 
the  first  multivariable  control  study  of  the  column.  The 
model  was  derived  from  non-steady  state  material  and  energy 
balances  linearized  around  the  steady  state  operating  con¬ 
ditions.  The  model  was  put  in  the  state  space  form  with 
twenty  states  being  the  liquid  enthalpies  and  the  liquid 
flows  from  each  stage.  Since  the  control  law  of  interest 
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required  that  all  states  be  measured  it  was  found  necessary 
to  approximate  the  model  using  a  subset  of  the  states.  This 
could  have  been  accomplished  by  estimation  techniques  (21,23) 

49)  or  by  model  reduction  (1,26,52).  He  used  Anderson  s  (1) 
method  to  reduce  the  model  to  second  order  retaining  only 
the  terminal  compositions  that  were  being  measured  and 
controlled.  The  reduction  algorithm  discretized  the  model 
enabling  control  to  be  carried  out  using  a  digital  computer. 
McGinnis  used  a  control  scheme  that  contains  both 
feedforward  and  feedback  elements  and  does  not  require 
explicit  decoupling  as  interactions  are  taken  into  account 
in  the  calculations.  Although  initial  runs  gave  unsatis¬ 
factory  oscillatory  response  due  to  the  noise  level  in  the 
measurement  of  the  bottoms  composition,  McGinnis  intuitively 
adjusted  the  values  of  elements  in  the  control  law  to  obtain 
satisfactory  control.  Comparison  of  the  resulting  control 
with  that  achieved  by  Berry,  however,  showed  that  the  control 
behavior  obtained  in  the  previous  study  was  superior.  This 
was  attributed  to  the  fact  that  McGinnis  did  not  account  for 
time  delays  in  his  model. 

Other  control  schemes  have  been  developed  which  could 
be  used  to  control  the  distillation  column.  Optimal  control 
by  definition  should  give  the  best  control,  however  because 
of  its  dependence  on  the  model  and  high  gains  often  encountered 
optimal  control  may  not  the  best  choise.  Hu  and  Ramirez 

ented  control  laws  for  the  optimal  control  of 


(18)  have  pres 
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distillation  columns. 

Rosenbrock  (45)  proposed  a  method  of  transforming  the 
linear  state  space  model  into  a  non-interacting  set  of 
modes.  This  new  system  can  then  be  used  to  calculate  control 
without  the  need  for  decoupling.  The  dependence  on  having 
an  accurate  model  is  again  the  major  shortcoming  for  this 
approach.  The  fact  that  the  physical  variables  are  not 
being  dealt  with  directly  adds  a  problem  in  visualizing  the 
actual  system.  More  recent  treatments  of  modal  control  have 
been  presented  by  Gould  (12)  and  by  Ellis  and  White  (9,10,11). 

The  origin  of  the  control  technique  used  in  this  work, 
Greenfield  and  Ward’s  method  which  they  call  ’’structural 
analysis",  makes  use  of  the  "structural"  information  avail¬ 
able  in  the  linearized  material  and  energy  balances.  The 
structural  information  they  use  is  contained  in  the  state 
space  model  and  is  simply  the  relationships  between  the 
internal  variables  (or  states)  that  are  not  given  explicitly 
by  transfer  function  models.  The  method  is  based  on  the 
"Invariance  Principle"  (17,36)  and  gives  complete  decoupling 
and  feedforward  control  action  with  allowance  for  the 
addition  of  feedback  loops.  Bollinger  and  Lamb  (3)  proposed 
a  similar  control  scheme  based  on  a  transfer  function  model 
form.  The  results  of  these  two  methods  and  one  called 
V-form  analysis  proposed  by  Mesarovic  (29)  have  been  compared 
in  a  series  of  articles  by  Rajagopolan  and  Sheshadri  (37,38, 
39)  for  a  continuous  stirred  tank  reactor.  It  was  found  that 
when  the  chemical  reaction  in  the  reactor  was  second  order 
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the  resulting  control  law  from  structural  analysis  required 
simple  gains  while  the  two  transfer  function  techniques 
required  lead-lag  controllers.  The  simplicity  of  the 
controllers  obtained  by  structural  analysis  is  due  to  the 
extra  relationships  between  the  states  that  are  available. 

It  is  on  the  basis  of  the  simpler  control  law  that  the  extra 
effort  required  to  get  the  linear  state  space  model  can  be 
justified . 

The  removal  of  the  interactive  effects  by  Greenfield 
and  Ward's  method  includes  both  feedforward  or  input  de¬ 
coupling  and  feedback  or  output  decoupling.  The  effective¬ 
ness  of  feedback  decoupling  has  been  discussed  by  Mitchell 
and  Webb  (30).  They  showed  that  adding  decoupling  does  not 
always  improve  control.  Niederlinski  (32)  presented  some 
guidelines  to  help  choose  between  decoupling  and  interactive 
control.  It  is  easier  to  implement  servo-control  on  de¬ 
coupled  loops.  Interactive  control  may  give  better  results 
for  regulatory  control  but  tuning  of  individual  loops  is 
difficult . 

The  final  feature  of  Greenfield  and  Ward's  method  is 
the  allowance  made  for  the  addition  of  feedback  loops.  The 
scheme  compensates  for  interaction  allowing  the  feedback 
loops  to  be  treated  as  a  single  variable  control  problem. 

The  advantage  of  this  feature  is  apparent  as  it  allows  on¬ 
line  tuning  which  permits  the  designer  to  compensate  for 
model  inadequacy. 
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4.2  Mathematical  Development  of  the  Control  Scheme 

The  response  of  a  process  such  as  distillation  can  be 
represented  by  a  set  of  first  order  non-linear  differential 
equations.  These  equations  can  be  linearized  by  retaining 
linear  terms  of  a  Taylor  series  expansion  and  the  resulting 
set  of  equations  can  be  put  in  the  state  space  form.  Using 
such  a  set  of  equations  Greenfield  and  Ward  have  developed  a 
scheme  to  decouple  the  system  and  introduce  feedforward 
control.  The  motivation  for  using  the  state  space  model 
rather  than  the  transfer  function  form  lies  in  what  Green¬ 
field  and  Ward  call  the  "structural"  information.  This  is 
the  equations  giving  the  response  of  the  internal  "structural" 
variables  or  simply  the  states.  These  relationships  simplify 
the  derivation  of  an  invariant  control  law.  The  "Invariance 
Principle"  as  discussed  by  Petrov  (36)  and  by  Haskins  and 
Sliepcevich  (17)  is  the  principle  behind  Greenfield  and 
Ward’s  control  scheme. 

The  development  shown  here  includes  first  the  continuous 
system  formulation  of  the  control  law  for  the  case  where  some 
states  are  not  measured  and  for  the  case  where  alH  states  are 
measured  or  independently  estimated.  The  discrete  control 
law  is  formulated  for  the  case  where  all  states  are  measured 
since  it  is  necessary  that  all  states  be  available. 

4.2.1  Continuous  Formulation 

The  linear  state  space  model  can  be  expressed  as 

=  Az  +  Bu  +  Cd 


z 


4.1 
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where  z_ 

u 

d 

A ,  B ,  C 


a  vector  of  the  states 
a  vector  of  the  control  inputs 
a  vector  of  the  disturbance  inputs 
model  matrix,  control  matrix  and 
disturbance  matrix  which  are  the 
coefficient  matrices. 


This  model  is  in  linear  perturbation  form  so  taking 
the  Laplace  transform  gives 


where 


s  Z  =  A  Z  +  BU  +  CD  4.2 

Z,U,D  =  the  Laplace  transforms  of 
£,u  and  d  respectively 
s  =  the  independent  variable  (Laplace) 


The  states  can  be  divided  into: 

Z^  -  controlled  states  (measured  or  estimated) 
Z2  -  measured  or  estimated  states  not 
controlled 

Z g  -  unmeasured,  uncontrolled  and  unestimated 
states 

By  partitioning  the  model,  control  and  disturbance  matrices, 
equation  4.2  can  be  written  as: 
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Note  that  it  is  assumed  that  the  controlled  states  Z_-^  and 
the.  manipulated  variables  U  are  ordered  such  that  the  first 
manipulated  input  will  be  used  for  feedback  control  of  the 
first  state  and  the  second  for  the  second  and  so  on  to  the 
last  controlled  state.  Also,  in  this  model  it  is  necessary 
that  the  number  of  independent  controlled  outputs  or  states 
be  less  than  or  equal  to  the  number  of  independent  mani¬ 
pulated  inputs  (13)  or: 


dim  (Z,^)  <  dim  (U) 


where  dim  (Z^)  -  the  dimension  of  the 


vector  Z 


-1 


In  the  case  where  the  system  does  not  facilitate  direct 
measurement  of  the  states,  it  is  generally  possible  to 
express  the  measured  outputs  as  a  linear  combination  of  the 
states  or: 


y 


Z 

[Hi  H2] 


4.4 


where  y  -  output  vector 

—l5 —2  “  output  matrices. 


In  this  relationship  it  is  necessary  that 


dim  (y)  =  dim  (Z^)  +  dim  (Z^) 


so  that  the  matrix  [H^  H^]  is  square  and: 


[Hi  H2]  1  x 


4.5 
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The  relationships  given  in  equations  4.3  and  4.4  are 
represented  in  the  signal  flow  graph  given  in  Figure  4.1. 

In  addition  to  the  dimensional  restrictions  already 
stated,  two  further  requirements  must  be  met  in  order  to 
fulfil  the  invariance  criterion: 


1.  Manipulated  inputs  are  assumed  unbounded 
under  control  demands. 

2.  No  pure  time  delays  are  present  in  the 
system. 

The  first  step  in  development  of  the  control  algorithm 
is  the  elimination  of  from  the  expression  for  s  Zq .  It 
follows  from  equation  4.3  that: 

s  2q  =  Aq  Zq  +  A2  Z2  +  A3  Z3  +  Bx  U  +  Cx  D  4.6 

In  order  to  eliminate  Z^  it  is  necessary  to  solve  for  Z^ 
from  the  equation: 


s  Z3  =  Ay  Zq  +  Ag  Z2  +  Ag  Z3  +  B3  U  +  C3  D  4.7 

Solving  for  Z^  gives 

Z3  =  (Si-  Ag)-1  (A7  Zx  +  Ag  Z2  +  Bg  U  +  C3  D)  4.8 


where  I  -  identity  matrix  of  appropriate  dimension 


Let 


G  =  (s  I  -  Ag) 


and  substitute  equation  4.8  into  equation  4.6. 
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FIGURE  4.1:  Signal  Flow  Graph  of  the  Model  Equations. 
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s  Zp  -  (Ap  ^3  ^  Ay)Z-^  +  ( Ag  +  Ag  _G  Ag )  Zg 
+  (ll  +  A3  G  |3)U  +  (C1  +  Ag  G  Cg)D 


Let 


4.9 


L  —  Ap  +  Ag  G  Ay 

M  =  Ag  +  Ag  G  Ag 

P  =  *!  +  Ag  G  Bg 

S  =  £1  +  i  £3 


4.10 


Separate  matrix  L  into  two  matrices  such  that  one  matrix, 
,  is  diagonal  containing  the  main  diagonal  of  L  and  the 
other,  L  ,  consists  of  the  original  matrix  with  the  diagonal 
set  to  zero  such  that 


L  ,  +  L 
=d  =p 


L 


4.11 


Substituting  relations  4.10  and  4.11  into  4.9  gives 


s  Z,  =  (L,  +  L  )Z,  +MZo+PU+0D 
—1  =d  =p  —1  =  —2  =  —  —  — 


4.12 


Rearranging  gives 


-1 


Zp  =  (s  I  -  Ld)  (L  Zp  +  M  Z2  +  P  U  +  £  D)  4.13 


This  represents  a  set  of  linear  equations  of  the  form 


Zi  =  fp  (Zp,  Z2,  U,  D) 


4.14 


Z.  -  the  ith  element  of  Zn 
1  —1 

Zp  -  vector  Zp  with  element  i  deleted 


where 


< 
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For  invariance  =  0  in  equation  4.14  or  equation  4.13 


becomes 

0  - 

:  (s  I  -  Ld)  1  CLp  Z±  +  M  Z2  +  P  U  +  g  D)  4.15 

where  0  - 

-  null  matrix 

which  will  be  satisfied  if 


Separate 

LpZ1+MZ2+PU+QD=0  4.16 

P  into  two  matrices  as  was  done  with  L  to  get 

PJ  +  P  =  P  4.17 

=d  =p 

and  substituting  into  equation  4.16  gives 

ip  iq  +  M  Z2  +  Pd  U  +  Pp  U  +  £  D  =  0  4.18 

which  upon  rearrangement  yields 


U  =  -P"1  (L  Z-.  +  M  Z0  +  P  U  +  a  D)  4.19 

-  =d  =p  -1  =  -2  =p  -  -  - 

which  is 

of  the  form 

Up  =  f±  (Z”,  Z2,  u"5  D)  4.19a 

where 

U.  -  element  of  U 

l  — 

U  -  vector  U  with  element  i  deleted 

Equation  4.19  gives  the  relation  used  for  feedforward  plus 
decoupling  control  as  obtained  by  Greenfield  and  Ward. 
Substituting  equation  4.5  gives 


U  =-£d1  ([L  ^][H1  u2l  1  z  +  £p  U  +  2  D) 


4.20 
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The  control  matrices  thus  are: 


where 


II* 

* 

* 

=  -Zl1  2, 

£dc 

--  -fd1[iP 

M][Hi  H2]  1 

4.21 

£uc 

=  -P^P 

=d  =p 

=FF 9 

£dc’  =uc 

control  matrices  - 

feedforward , 

output  decoupling,  input  de¬ 
coupling. 


The  control  system  is  represented  with  a  signal  flow  graph 
in  Figure  4.2. 

At  this  point  feedback  control  action  can  be  included 
in  the  control  function: 


Let  V  -  the  feedback  portion  added  V.  -  f.  (Z.) 

—  ^  111 

T 

(for  P  +  I  feedback  V.  =  K  Z .  +  KT  f  Z.dt  at  t=T) 

i  p  i  I  ^  i 


t=  0 


thus  adding  V  to  equation  4.20  gives 


U  =  -P,1  (L  Z,  +  M  Z0  +  P  U  +  0  D)  +  V 

—  =d  =p  —1  =  —2  =p  —  ^  —  — 


4.22 


For  the  case  where  all  states  are  measured  equation 
4.2  can  be  written  as 
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where  the  states  are  as  defined  previously. 
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FIGURE  4.2:  Signal  Flow  Graph  of  the  System  with 

Feedforward  and  Decoupling  Control. 
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It  thus  follows  that  s  Z^  is 

s  Zx  =  A1  Zx  +  A2  Z2  +  J1  U  +  C]_  D  4.24 

Forming  the  diagonal  matrix  and  the  matrix  consisting 
of  A^  with  the  main  diagonal  set  to  zero  such  that 


4.25 


Substituting  into  equation  4.24  and  rearranging  for  gives 
Zx  =  (si-  Ad)_1  (A  Zd  +  A2  Z2  +  B1  U  +  C±  D)  4.26 


As  before  invariance  holds  if 


Ap  Zd+A2  Z2+  U  +  C1  D  =  0  4.27 

Separation  of  into  a  diagonal  matrix  and  a  matrix 
with  a  zero  main  diagonal  such  that 


B  +  B,  =  Bn  4.28 

=p  =d  =1 


and  substitution  into  equation  4.27  and  rearranging  yields: 


U 


(£p  h.  + 


a2  z2 


+  B  U  +  C-,  D) 
=P  “  =1  “ 


4.29 


This  is  the  control  function  in  the  form  presented  by  Green¬ 
field  and  Ward.  Unlike  the  previous  case  where  some  states 
were  eliminated  algebraically  this  result  has  constant  control 
matrices.  An  equivalent  expression  for  U  is  obtained  by 


solving  equation  4.27 


U 


(A>  h. 


+  C±  D) 


4.30 


. 
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Expressing  equation  4.30  in  terms  of  the  output  gives: 


U  =  -B11  ([A  A2][H1  H2]  1  y  +  C1  D) 


4.31 


giving  control  equations 


-B-,1  [=p  =2]C£l  =2  ^ 


4.31a 


The  difference  between  these  results  is  that  equation  4.30 
unlike  equation  4.29  does  not  contain  explicit  decoupling  of 
the  control  inputs.  This  becomes  important  when  the  feed¬ 
back  portion  of  the  control  is  added.  Adding  feedback  control 
to  equation  4.29  gives: 


U  =  -B,1  (A_  Z-.  +  A0  Z0  +  B  U  +  C-,  D)  +  V 
—  =d  =p  —1  =2  —2  =p  —  =1  —  — 


4.32 


while  inclusion  with  equation  4.30  yields 


U  =  -B11  (A  Zx  +  A2  Z2  +  C1  D)  +  V 


4.33 


To  illustrate  the  difference  a  third  order  example  is 
presented  in  Appendix  C.  The  control  function  as  given  in 
equation  4.31  is  illustrated  in  Figure  4.3.  The  advantage 
of  this  form  is  that  when  control  is  carried  out  by  a  digital 
computer  it  is  unnecessary  to  solve  the  set  of  equations  for 
each  control  output  by  numerical  means  as  would  be  the  case 
using  equation  4.29  Cor  4.32)  because  U  is  given  explicitly. 
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FIGURE  4.3:  Signal  Flow  Graph  of  System  Incorporating 

Simplified  Control  Scheme. 
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4.2.2  Discrete  Formulation 

The  previous  formulation  of  the  control  function  can 
also  be  adapted  for  the  discrete  version  of  the  linear  state 
space  model  although  an  additional  restriction  must  be  met. 

In  the  discrete  formulation  the  elimination  of  unmeasured 
states  from  the  model  cannot  be  done  directly  so  that  states 
that  cannot  be  measured  or  estimated  must  be  removed  by  some 
method  such  as  mod  3.1  reduction. 

The  discrete  model  equivalent  to  the  continuous  linear 
state  space  model  used  in  the  previous  section  can  be  ex¬ 
pressed  as 

zCOc+DT]  =  (T)  zCkT)  +  £(T)  u(kT)  +  ACT)  d(kT)  4.34 

where  T  -  the  discrete  time  interval 

k  -  the  integer  number  of  intervals 

z(kT)  -  vector  of  states  at  time  kT 

u(kT)  -  vector  of  control  inputs  at  time  kT 

d(kT)  -  vector  of  disturbance  inputs  at  time  kT 

^(T)  -  fundamental  matrix 

j|(T)  -  control  matrix 

A(T)  -  disturbance  matrix 

For  simplicity  the  symbol  for  the  discrete  time  interval  T 
will  be  omitted  from  the  following  derivation. 

Dividing  the  states  into  those  controlled  z_-^  and  those 
not  controlled  allows  the  model  to  be  given  as : 
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zx(k+l) 

±i 

±2 

z^(k) 

+ 

u(k)  + 

'  H 

<111 

1 _ 

^2  Ck+1) 

is 

i4 

z2  (k) 

—2 

1 

ll> 

K> 

1 _ 

The  relation  governing  the  response  of  z^  is  given  by 

z^Ck+l)  =  z^Ck)  +  z2Ck)  +  u(k)  +  A1  d(k)  4.36 

The  first  step  in  finding  the  control  law  is  to  divide  , 
into  the  diagonal  matrix  and  the  matrix  J>p  with  the 
original  elements  except  for  the  main  diagonal  which  is  set 
to  zero,  such  that 


ip  +  id  =  ii  4-37 

Substituting  into  equation  4.36  gives 

z^(k+l)  =  (_£  +  z^(k)  +  ^2  z^Ck)  +  iq  H.^) 

+  Ax  d (k)  4.38 

which  upon  rearrangement  yields 

z^Ck+l)  -  z-^Ck)  =  z_^(k)  +  ^2  z_2(k)  +  u(k) 

+  d( k)  4.39 

This  equation  has  the  form 

z.(k+l)  -  d>,  z.(k)  =  f.  (z_n  ,  z?(k),  uCk),  dCk))  4.40 

1  Cl  •  •  1  11^ 

11 

where  z.  -  ith  element  of  zn 

l  —1 

Tt 

z^  -  vector  z^  with  element  i  deleted 

For  invariance  at  the  sampling  instances  in  the  control 
variable  z^  it  is  necessary  that 
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z^Ck)  =  Zq(k+1)  =  z(k+n)  =  0,  for  n  =  any 
positive  integer. 

This  means  that  the  left  hand  side  of  equation  4.40  is  zero 
so  equation  4.39  becomes: 


£  =  z^Ck)  +  +  J=i  +  Ai  £(k) 


4.41 


D 


ividing  into  the  diagonal  matrix  and  the  matrix  with 


zero  on  the  main  diagonal  j3^  such  that 


Ip  +  Id  "  II  4 

and  substituting  this  into  equation  4.41  and  rearranging 
gives 


uCk)  =  -Cj$d)  1  £q(k)  +  £2^)  +  Ip 

+  Aq  d(k)) 


4.43 


The  control  matrices  can  be  given  as 

A 

Iff  =  "Id  Ii 

Sue  "  "Id  Ip 

=DC  =  “Id  [ipl2][Sl  =2  ^ 

AAA 

where  Iff5  — UC 5  -DC  ~  discrete  control  matrices  for 

feedforward,  input  and  output 
decoupling  respectively. 


This  is  the  control  function  for  the  discrete  case  and  the 
signal  flow  graph  for  this  is  given  in  Figure  4.4  with 
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output  feedback  where 

y(k)  =  [H^  H2] 


Ck) 


so 


^Ck) 

z_2 


z_2  (k) 


[Hi  H,]’1  yCk) 


where  dim  y(k)  =  dim  Cz^Ck)) 


4.44 


4.45 


As  in  the  continuous  case  the  analytic  solution  of 
equation  4.41  gives  an  approximate  solution  lacking  explicit 
control  input  decoupling.  This  is  given  as 

uCk)  =  -j}”1  z^Ck)  +  <^2  220<)  +  A1  dCk))  4.46 

The  motivation  for  using  this  form  is  greater  for  the 

A 

discrete  case  as  the  algebraic  loop  including  must  be 

solved  by  an  iteritive  method  at  each  sample  interval.  The 
addition  of  feedback  control  such  as  the  discrete  approxi¬ 
mation  of  P  +  I  control  given  by 

k 

W . (k)  =  K  z.(k)  +  Kt  Z  (z.(j)  T)  4.47 

1  P  1  1  j=o  1 

For  the  case  with  explicit  decoupling 

u(k)  =  jl”1  l$_  z,  (k)  +  j>0  z0(k)  +  8  u(k) 

a  p  -i  z  z  P  4.48 

+  d Ck) ]  +  W(k) 
and  for  the  simplified  control  as 

uCk)  -  z^(k)  +  z_2^^  +  =q  d)  +  W(k)  4.49 

This  is  illustrated  with  a  signal  flow  graph  in  Figure  4.5. 
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FIGURE  4.4:  Signal  Flow  Graph  of  Discrete  System  with 

Feedforward  and  Decoupling  Control. 


FIGURE  4.5: 


Signal  Flow  Graph  of  Discrete  System  with 
Simplified  Feedforward  and  Decoupling 
Scheme  with  Feedback  Loop  Shown. 
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4.3  Implementation  of  the  Control  Scheme 

The  method  of  deriving  feedforward  and  decoupling 
control  for  multivariable  systems  was  applied  to  the  pilot 
scale  distillation  column  at  the  University  of  Alberta.  The 
design  procedure  used  to  implement  control  included  first 
the  calculation  of  the  control  law  from  the  model  of  the 
system  then  a  simulation  study  followed  by  verification  of 
the  effectiveness  of  the  scheme  on  the  experimental  equipment. 

4.3.1  Simulation 

The  model  derived  in  the  previous  chapter  was  used  to 
calculate  the  control  laws.  The  linear  continuous  version 
of  the  distillation  column  model  was  used  with  equation  4.33 
to  obtain  the  corresponding  control  matrices.  It  was  found 
that  the  control  matrices  were  singular  when  the  control 
variables  were  the  terminal  compositions  or  the  first  and 
tenth  states.  This  was  due  to  a  row  of  zeros  in  ^  in 
equation  4.33  corresponding  to  the  lack  of  direct  effect  of 
the  inputs  on  the  overhead  composition.  From  the  model 
matrices  given  in  Appendix  B  it  can  be  seen  that  the  only 
variable  that  affects  the  overhead  composition  is  the  com¬ 
position  on  the  adjacent  stage.  It  is  therefore  logical  to 
use  this  as  the  control  variable  in  place  of  the  overhead 
composition. 

The  control  matrices  resulting  from  the  calculation 
using  the  bottoms  composition  and  the  composition  of  the 
liquid  on  stage  nine  (eighth  tray)  are  given  in  Appendix  D. 
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Due  to  the  tridiagonal  form  of  the  model  of  the  system  it 
is  necessary  to  measure  only  states  two,  eight  and  ten  in 
order  to  decouple  the  system. 

The  response  of  the  system  with  the  calculated  control 
to  a  ten  percent  increase  in  feed  was  predicted  using  CSMP 
III  (56).  This  is  an  IBM  program  developed  for  solving 
differential  equations.  The  predicted  response  using  the 
linear  tenth  order  continuous  model  gave  almost  no  deviation 
from  steady  state  as  can  be  seen  from  the  final  deviation 
given  in  Table  4.1  run  4.1.1.  The  control  function  used 
here  does  not  include  feedback  but  from  the  results  they  do 
not  appear  to  be  required.  In  order  to  give  a  more  realistic 
prediction  the  control  matrices  were  used  with  the  non-linear 
model  of  the  system.  Even  though  the  results  were  not  as 
good  as  the  previous  ones  they  still  indicated  that  almost 
invariant  control  of  the  terminal  compositions  can  be 
expected  as  seen  by  the  final  deviation  given  in  Table  4.1 
run  4.1.2. 

Since  the  actual  control  function  for  the  physical 
system  is  carried  out  at  discrete  one  minute  intervals  by  a 
program  on  the  IBM  1800  digital  computer  the  discrete  control 
law  was  found.  In  order  to  do  this  it  is  necessary  to  dis¬ 
cretize  the  system  model.  The  discrete  model  was  calculated 
using  a  program  in  GEMSCOPE.  GEMSCOPE  is  a  collection  of 
programs  written  in  the  Department  of  Chemical  Engineering 
at  the  University  of  Alberta  to  carry  out  control  problem 
calculations.  These  are  described  by  Wilson  (51)  who  linked 
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TABLE  4.1 

DEVIATION  OF  TERMINAL  COMPOSITIONS  PREDICTED  FROM  TENTH  ORDER  LINEAR  CONCENTRATION 
MODEL  WITH  FEEDFORWARD  DECOUPLING  FUNCTION  FOR  A  TEN  PERCENT  STEP  INCREASE  IN  FEED  FLOW 
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them  into  a  convenient  form  for  the  IBM  360  operating  under 
the  MTS  system.  The  resulting  tenth  order  model  matrices  of 
the  discrete  linear  model  with  a  one  minute  sample  time  are 
given  in  Appendix  D.  The  response  of  this  model  was  compared 
to  that  of  the  continuous  version  for  a  ten  percent  increase 
in  feed  and  found  to  agree  well.  The  control  matrices  were 
calculated  from  the  discrete  model  according  to  equation 
4.46  and  the  response  with  the  control  to  a  ten  percent 
increase  in  feed  predicted  using  a  program  in  GEMSCOPE. 

The  control  matrices  are  given  in  Appendix  D  and  the  final 
steady  state  of  the  control  run  in  Table  4.1,  run  4.1.3. 

The  discrete  model  as  used  to  calculate  the  controlled 
response  resulted  in  almost  invariant  response  and  so  this 
form  of  control  would  seem  best.  However,  in  discretizing 
the  model  the  tridiagonal  form  of  the  model  matrix  was  lost 
and  so  the  decoupling  portion  required  that  all  the  states 
be  measured.  This  was  thought  wasteful  of  computer  time  and 
storage  space. 

In  order  to  determine  if  simplified  control  was 
feasible  it  was  decided  to  reduce  the  model.  The  GEMSCOPE 
program  contains  two  different  model  reduction  routines, 
modal  reduction  proposed  by  Marshall  (26)  and  least  squares 
reduction  proposed  by  Anderson  (1)  and  modified  by  Wilson 
(52).  The  tenth  order  model  was  first  reduced  to  a  fourth 
order  model  retaining  the  terminal  compositions  and  the 
liquid  compositions  of  the  second  and  ninth  stages.  The 
results  from  both  reduced  models  (one  from  least  squares  and 
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one  by  modal  reduction)  gave  excellent  final  steady  state 
agreement  with  the  predicted  final  steady  state  of  the  tenth 
order  discrete  model  as  can  be  seen  in  Table  4.2,  runs  4.2.2, 
4.2.3  and  4.2.4.  However,  the  dynamic  agreement  of  the 
model  reduced  by  the  modal  technique  was  by  far  the  better 
of  the  two  as  can  be  seen  in  Figures  4.6  and  4.7.  Reduction 
to  second  order  gave  results  similar  to  those  of  the  fourth 
order  models  as  shown  in  Table  4.2  and  Figure  4.7. 

The  control  functions  were  calculated  using  the  fourth 
and  second  order  reduced  models  using  equation  4.46  and  the 
resulting  control  matrices  are  presented  in  Appendix  D. 

The  predicted  controlled  response  to  a  ten  percent 
increase  in  feed  from  the  tenth  order  model  was  calculated 
using  the  fourth  and  second  order  control  matrices.  As  can 
be  seen  from  the  final  steady  states  given  in  Table  4.1, 
runs  4.1.4  and  4.1.5,  the  terminal  compositions  did  not 
deviate  significantly  indicating  that  essentially  invariant 
control  was  achieved. 

It  was  noted  that  the  feedback  decoupling  matrix  had 
a  negative  element  corresponding  to  feedback  from  the  over¬ 
head  composition  to  the  overhead  flowrate  in  the  fourth 
order  control  law.  For  single  loop  control  this  would 
result  in  unstable  operation.  A  simulation  run  was  carried 
out  with  this  element  set  to  zero.  This  resulted  in  almost 
no  change  in  the  predicted  final  steady  state  as  can  be  seen 
in  Table  4.1,  run  4.1.5. 


. 


PREDICTED  OPEN  LOOP  DEVIATION  OF  TERMINAL  COMPOSITIONS  TO  A  TEN  PERCENT 

INCREASE  IN  FEED  FLOW  USING  REDUCED  MODELS 
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FIGURE  4.6: 


Experimental  Open  Loop  Response  to  a  Ten 
Percent  Step  in  Feed  (2.4-2.64  LB/MIN) 
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FIGURE  4.7:  Comparison  of  Predicted  Response  of  Tenth 
Order  Linear  Discrete  Model  and  Second  and 
Fourth  Order  Reduced  Model  for  a  Ten  Percent 
Step  in  Feed 
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Consideration  of  the  results  from  the  discrete  sim¬ 
ulation  reveals  that  the  predicted  responses  of  the  tenth 
order  model  with  the  various  control  matrices  exhibited 
final  offsets  which  are  far  smaller  than  can  be  measured 
experimentally.  Thus  it  does  not  appear  to  matter  which  set 
of  control  matrices  is  chosen  from  a  control  quality  stand¬ 
point.  However,  two  other  features  must  be  considered  in 
implementing  the  scheme  on  the  actual  system.  First, 
measurement  noise  and  unreliability  present  in  the  physical 
system  can  render  a  sensitive  system  such  as  one  with  large 
feedback  elements  inoperative.  Secondly,  the  model  of  the 
process  is  not  perfect,  for  example,  the  distillation  column 
is  quite  non-linear  while  the  control  law  formulation  is 
based  on  a  linear  system  representation  which  means  offset 
can  be  expected  on  any  disturbance. 

In  the  simulation  runs  the  addition  of  the  feedback 
loops  was  not  needed  because  the  feed  disturbance  introduced 
could  be  adequately  attenuated  by  the  feedforward  action  of 
the  control  law.  However,  as  feedback  will  probably  be 
required  in  the  experimental  study  some  simulation  was 
carried  out  with  integral  feedback  loops  to  determine  what 
the  effect  will  be.  In  order  to  study  the  effect  of  feed¬ 
back  the  feedforward  control  was  removed  and  integral  feed¬ 
back  loops  were  added  to  the  overhead  and  bottoms  composi¬ 
tions.  The  results  of  the  tests  made  are  shown  in  Table  4.3. 
Runs  4.3.1  and  4.3.2  show  the  effect  of  removing  the  feed¬ 
forward  action  leaving  only  the  decoupling  function.  The 
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TABLE  4.3 

PREDICTED  DEVIATIONS  OF  TERMINAL  COMPOSITIONS  FOR  A  TEN 
PERCENT  INCREASE  IN  FEED  FLOW  USING  THE  TENTH  ORDER 
CONCENTRATION  MODEL  WITH  DECOUPLED  FEEDBACK  CONTROL 


Run 

Deviation  of 

Bottoms 

Composition  at 

Kl 

100  Min  (Wt.  Pet. 

Overhead 

MEOH) 

KI 

4.3.1* 

.  0029 

0 

-.00303 

0 

4.3.2 

.  0032 

0 

-.00056 

0 

4.3.3 

.  0023 

100. 

.  0002 

1. 

4.3.4 

.  0033 

1. 

.5  x  10"5 

1. 

4.3.5 

.  0034 

0 

.3  x  10"5 

1. 

4.3.6 

.000067 

1000  . 

.9  x  10"4 

1. 

4.3.7 

.25  x  10”5 

10000. 

-4 

.2  x  10 

1. 

* 


Negative  element  in  feedback  decoupling  matrix  retained 
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resulting  final  steady  state  is  different  from  the  open 
loop  run  for  which  the  final  steady  state  is  given  in 
Table  4.2,  run  4.2.1.  Removal  of  the  negative  element 
in  the  decoupling  matrix  corresponding  to  feedback 
from  the  distillate  composition  to  the  distillate  flowrate 
decreases  the  final  steady  state  offset  of  the  distillate 
composition,  as  can  be  seen  from  Table  4.3,  and  furthermore, 
the  overall  control  performance  is  improved.  It  is  interest¬ 
ing  to  note  that  both  increasing  and  decreasing  the  bottoms 
integral  constant  from  100  resulted  in  a  decrease  in  the 
final  overhead  composition  offset.  The  bottoms  composition 
deviation  at  final  steady  state  decreases  up  to  an  integral 
constant  of  ten  thousand  at  which  point  some  evidence  of 
oscillation  was  present  but  the  system  remained  stable. 

The  noise  level  in  the  measurement  devices  is  an 
important  feature  of  the  experimental  equipment.  The  gas 
chromatograph  used  to  measure  bottoms  composition  has  a  very 
high  noise  level  as  can  be  seen  from  the  steady  state  data 
plotted  in  Appendix  A.  There  is  a  time  delay  of  three  minutes 
for  the  analysis  in  addition  to  the  delay  resulting  from 
exponential  filtering.  Although  the  tray  temperatures 
which  are  used  to  estimate  liquid  compositions  also  have 
relatively  high  noise  levels,  they  give  better  results  than 
the  gas  chromatograph. 

The  control  law  derived  from  the  fourth  order  discrete 
model  was  chosen  for  the  experimental  study.  It  has  the 
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advantage  of  allowing  some  internal  state  measurements  to 
supplement  the  terminal  composition  measurements  while  still 
having  small  dimension  thus  saving  computer  time  and  storage. 

The  success  of  the  simulation  study  still  leaves 
unanswered  the  question  of  performance  of  the  continuous 
control  function  with  the  experimental  equipment  and  also 
what  will  be  the  effect  of  system  noise  upon  performance. 
Consequently,  the  experimental  study  was  undertaken  to 
attempt  to  provide  answers  to  these  questions. 

4.3.2  Experimental 

The  object  of  the  control  study  was  to  obtain  a  simple 
efficient  control  scheme  for  the  pilot  plant  distillation 
column. 

The  first  phase  of  the  experimental  study  involved 
testing  the  continuous  control  law.  The  system  was  con¬ 
figurated  as  shown  in  Figure  4.8.  On  initiation  of  control 
at  steady  state  monitoring  the  terminal  compositions  revealed 
oscillations  of  increasing  magnitude  to  the  system  limits 
before  the  feed  disturbance  was  introduced.  The  system  was 
unstable  so  it  was  necessary  to  put  the  control  law  in  the 
discrete  form. 

Control  of  the  column  was  then  attempted  using  the 
simplest  discrete  control  law  resulting  from  the  use  of  the 
second  order  model  representation.  It  was  found  that  the 
system  oscillated  due  to  the  high  noise  level  in  the  gas 
chromatograph  measurement  of  bottoms  composition.  The 
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FIGURE  4.8: 


Schematic  of  Column  Control  Configuration 
for  Continuous  Control  Law 
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system  was  stable  but  control  was  found  to  be  unsatisfactory. 

The  fourth  order  model  which  was  considered  the 
logical choise  was  used  to  find  the  control  law,  the  numeri¬ 
cal  values  for  which  are  given  in  Appendix  D.  A  simplified 
schematic  of  the  system  layout  is  shown  in  Figure  4.9.  The 
final  steady  state  deviations  of  the  controlled  variables 
are  given  in  Table  4.4  and  the  corresponding  feedback  con¬ 
stants  and  filter  constants  used  for  the  bottoms  composition 
measurement  are  given  in  Table  4.5. 

The  response  of  the  system  with  the  fourth  order  feed¬ 
forward  control  function  to  an  increase  and  decrease  of  ten 
percent  are  shown  in  Figures  4.10  and  4.11  respectively. 

The  disturbance  is  well  attenuated  but  a  final  offset  occurs 
in  both  runs . 

The  addition  of  feedback  control  with  decoupling 
action  resulted  in  the  improved  response  to  feed  flow 
changes  as  shown  in  Figures  4.12  and  4.13.  The  deviation 
in  both  compositions  appears  to  be  less  than  the  noise  level 
and  so  control  was  considered  to  be  satisfactory. 

Figures  4.14  and  4.15  show  the  effect  of  increasing 
the  integral  feedback  constants,  using  the  values  given  in 
Table  4.5,  where  it  appears  that  the  increase  in  integral 
action  gives  some  improvement  in  the  control.  To  test  this 
conclusion,  tests  were  made  using  a  series  of  random  dis¬ 
turbances.  The  results  as  shown  in  Figure  4.16  and  4.17 
reveal  that  the  higher  integral  constants  tend  to  make  the 
compositions  oscillate  and  for  these  larger  disturbances 
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PROCESS  LINES 
INSTRUMENT  LINES 


FIGURE  4.9:  Simplified  Schematic  of  Column  Control  for 

Discrete  Control  Law 


SUMMARY  OF  EXPERIMENTAL  RESULTS  FOR  CONTROL  TESTS 
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TABLE  4.5 

FILTERING  AND  FEEDBACK  CONTROLLER  CONSTANTS 
USED  IN  EXPERIMENTAL  TESTS 


Run 

Exponential  Filter 
Constant  + 

3Actual(min)  gDDC 

Integral  Controller  Feedback 
Constants  (sec) 

Bottoms  Loop  Overhead  Loop 

1 

0 

0 

0 

0 

2 

0 

0 

0 

0 

4 

0.5 

0.84315 

0 

0 

5 

0.5 

0.84315 

0 

0 

6A 

0.5 

0.84315 

0 

0 

7 

0.5 

0.84315 

100. 

1. 

8 

0.5 

0.84315 

100. 

1. 

9 

0.5 

0.84315 

100. 

1. 

11 

0.5 

0.84315 

100. 

1. 

13 

0.5 

0.84315 

100. 

1. 

14 

0.5 

0.84315 

1000. 

2. 

16 

.75 

0.92625 

1000. 

2. 

20 

.75 

0.92625 

1000. 

2. 

21 

.75 

0.92625 

100. 

1. 

+  Constant  used  in  bottoms  composition,  Xg ,  filter 
equation  xB(i)=xB(i-l)  (1-g)  +  ^(i)measured 


* 


^Actual 

eDDC 


-  Constant  based  on  one  minute  sample  time 

Constant  based  on  actual  sample  time  of 
16  seconds 


.1 

\ 

FLOW  rate:  CLB/MIN)  composition  CHANGE  (WT  PC) 

3*0  -1  0  1 
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FIGURE  4.10:  Experimental  Response  with  Decoupled  Feedforward 

Control  for  a  10%  Step  in  Feed  (2.4-2.64  LB/MIN) . 


FLOW  RATE  (LB/M IN)  CGMTOBITICN  CHANGE  (WT  PC) 

3-0  -1  0  i 
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FIGURE 


TIME-  in  minutes 

4.11:  Experimental  Response  with  Decoupled  Feedforward 

Control  for  a  -10%  Step  in  Feed  Flow 
(2.4-2.16  LB/MIN) . 


FLCW  RATE  (LB/MINI  COMPOSITION  CHANGE  CWT  PC) 

0.5  3*0  -101 


100 


Experimental  Response  with  Decoupled  Feed¬ 
forward  and  Feedback  Control  for  a  10% 

Step  in  Feed  (2.4-2.64)  ( Kj  =  [l  ,1 00] ) . 


FIGURE  4.1'2: 


FLOW  RATE  (LB/MIN)  COMPOSITION  CHANGE  (WT  PC) 

3.0  -1  0  1 
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Experimental  Response  with  Decoupled  Feed¬ 
forward  and  Feedback  Control  for  a  -10% 
Step  in  Feed  Flow  (2.4-2.16  LB/MIN) 
(£-£=[1,100])  . 


FIGURE  4.13: 


FLOW  RATE  (LB/M IN)  CQMPOSITICN  CHANGE  (WT  PC) 

3-0  -1  0  * 
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time  in  minutes 


FIGURE  4.14:  Experimental  Response  with  Decoupled  Feed¬ 
forward  and  Feedback  Control  for  a  10% 

Step  in  Feed  Flow  (2.4-2.64  LB /MIN) 

( K x  =  [ 2  ,100  0])  . 


FLOW  RATE  (LB/MIN)  COMPOSITION  CHANGE  <WT  PC) 

O.s  3*0  -10  1 
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Experimental  Response  with  Decoupled  Feed¬ 
forward  and  Feedback  Control  for  a  -10% 
Step  in  Feed  Flow  (2.4-2.16  LB/MIN) 

(K^.  =  [2,1000]) . 


FIGURE  4.15: 


FLOW  RATE-:  (LS-'MIN)  COMPOSITION  CHANGE'  CWT  PC) 

3*0  -10  1 
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o  cVi  js  yp.  i p.o 

timf:  in  minutes 


Experimental  Response  with  Decoupled  Feed¬ 
forward  and  Feedback  Control  to  a  Series 
of  Disturbances  ( K =  [1,100]  ) 


FIGURE  4.16: 


FLOW  RATE  (LB/MIN)  COMPOSITION  CHANGE  (WT  PC! 

3-0  -10  1 
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r 


of  Disturbances  ( K^= [ 2 , 100 D 3 ) 


FIGURE  4.17: 
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they  result  in  poorer  control. 

Figures  4.18  to  4.20  give  some  results  from  the  work 
of  Berry  and  McGinnis  for  comparison  purposes.  The  run  of 
Berry  is  comparable  to  results  obtained  in  this  study. 
McGinnis’  control  scheme  showed  oscillatory  behavior  due  to 
the  noise  in  the  bottoms  composition  measurement.  McGinnis 
used  a  second  order  discrete  control  law  which  was  comparable 
to  the  second  order  scheme  tried  in  this  study  but  did  some 
intuitive  ’’tuning”  to  make  it  work. 

Figure  4.21  shows  the  response  of  the  control  variables 
to  a  change  in  feed  enthalpy.  This  test  was  carried  out  to 
ensure  that  other  types  of  disturbances  would  be  compensated 
for  by  the  scheme  employed  in  this  study. 


I 

al  nx  e  :tnBrto  o 


COMPOSITION  CHANGE  (wt  %) 
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THIS  STUDY  (  F  =  2.40 — ►  2.88  Ib/min  ) 


BOTTOMS 


Experimental  Response  Obtained  by  Berry  With 
Decoupled  Feedback  Control  Compared  to 
Response  With  Decoupled  Feedforward  Feedback 
Control  for  an  Increase  in  Feed  Flow. 


FIGURE  4.18: 


FLOW  RATE  (LB/MIN)  COWFOSITICN  CHANGE  .  WT  PC) 

3-0  -1  0  .  f 
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c  •  • 

FIGURE  4.19:  Experimental  Response  Obtained  by  M  Ginnis 

with  Multivariable  Control  for  a  5%  Change 
in  Feed  Flow. 


FLOW  RATE  (LB/MIN)  CCMPOSITICN  CHANGE  CWT  PC) 

3*0-1  o  l 
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o  •  • 

FIGURE  4.20:  Experimental  response  Obtained  by  M  Ginnis 

with  Multivariable  Control  for  a  -5% 

Change  in  Feed  Flow. 


FLOW  RATE  (LB/M IN)  COMPOSITION  CHANGE  .  CWT  PC) 
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FIGURE  4.21:  Experimental  Response  with  Decoupled  Feed¬ 
forward  and  Feedback  Control  for  a  Change 
in  Feed  Temperature  from  161-141°F 
(^  =  [1,100]). 
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CHAPTER  5 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  following  conclusions  can  be  drawn  from  the  three 

phases  of  this  study. 

1.  A  direct  material  balance  control  system  utilizing 
distillate  flow  to  control  overhead  composition 
gives  better  response  when  the  distillation  column 

is  subjected  to  energy  changes  rather  than  flow  changes. 
This  should  be  true  for  all  disturbances  when  both 
terminal  compositions  are  under  control  with  the 
reboiler  duty  used  as  one  of  the  manipulated  variables 
for  composition  control. 

2.  For  changes  in  feed  flow,  an  energy  balance  control 
scheme  with  the  composition  controlled  by  the  reflux 
flow,  gives  superior  control  performance  to  the 
direct  material  balance  control  scheme. 

3.  The  response  of  the  distillation  column  was  predicted 
using  a  material  balance  model  in  the  linear  state 
space  form.  Although  the  actual  system  is  highly 
nonlinear,  for  control  purposes  the  linearized  model 
provides  an  adequate  prediction  even  though  the 
agreement  with  experimental  open  loop  responses  was 
not  very  good. 

4.  Feedforward  and  decoupling  control  matrices  as  proposed 
in  this  study  have  the  advantage  of  allowing  for 
feedback  control  tuning.  The  control  scheme  based  on 
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a  theoretical  model  can  be  implemented  with  on-line 
tuning  to  compensate  for  shortcomings  in  the  model. 

5.  By  eliminating  algebraic  loops  the  control  scheme  can 
be  implemented  with  a  digital  computer  and  gives  suf¬ 
ficient  decoupling  for  regulatory  action. 

6 .  The  use  of  a  reduced  model  to  find  the  control  scheme 
gave  a  reasonable  control  performance.  Restriction  of 
the  minimum  dimension  of  the  reduced  model  was  dictated 
by  the  measurement  noise  rather  than  system  dynamics . 

In  this  study  the  high  noise  level  in  the  bottoms 
measurement  made  it  necessary  to  retain  two  other 
states  in  addition  to  the  terminal  compositions. 

7.  Since  the  control  matrices  contain  only  proportional 
elements  the  continuous  version  of  the  control  law 
should  be  possible  to  implement  with  conventional 
pneumatic  or  electronic  instrumentation. 

The  following  are  some  recommendations  for  future 

work . 

1.  Time  delay  and  noise  associated  with  the  gas 

chromatograph  used  to  measure  the  bottom  composition 
are  a  major  source  of  disturbance  in  the  system. 
Improvement  in  the  equipment  is  an  area  that  requires 
some  work.  Possibilities  include  enlarging  sample 
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size  and  speeding  vaporization  of  the  sample  by 
increasing  temperature  or  gas  velocity  in  the 
sample  area. 

2.  In  addition  to  improving  the  measurement  equipment 

a  comparison  of  the  use  of  a  Kalman  filter  to  condition 
the  measurement  value  to  exponential  filtering  would 
be  relevant. 

3.  This  system  is  also  suitable  for  testing  the  application 
of  a  Smith  Predictor  to  lessen  the  effect  of  the 

time  delay  associated  with  sampling. 

4.  A  more  general  method  of  improving  control  would 

be  the  addition  af  a  state  estimation  technique  such 
as  the  steady  state  type  proposed  by  Weber  and 
Brosilow  (50).  This  could  be  compared  to  the 
dynamic  observer  of  Luenberger  (23).  Both  methods 
would  enable  use  of  a  tenth  order  concentration 
model  based  control  law  with  temperature  measurements 
used  to  estimate  the  states. 

5.  Model  reference  adaptive  control  because  of  its 

self  tuning  feature  is  suitable  for  use  with  a  simple 
model.  Also  the  continuous  tuning  feature  would 
provide  compensation  for  the  nonlinearity  of  distil¬ 
lation  columns.  The  pilot  scale  distillation  column 
is  suited  to  this  type  of  test. 
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NOMENCLATURE 


A 

ADC 


ARC 


B 

B 

C 

d 

D 

D 

DAC 


DDC 

F 

G 

h 

H 

H 

I 

I. 

3 

IAE 

IE 

K 


model  matrix 

analog  to  digital  converter  for  digital 
computer  input 

analyzer  recording  controller  generally 

capacitance  dynalog 

bottoms  flow  (lb/min) 

control  matrix 

disturbance  matrix 

T 

disturbance  vector  [F  Xp  hp] 

Laplace  transform  of  d 
distillate  flow  (lb/min) 

digital  to  analog  converter  for  digital 
computer  output 

Direct  Digital  Control  -  IBM  control  program 
feed  flow  (lb/min) 

(si  -  A)-1 

liquid  enthalpy  (BTU/lb) 

output  matrix 

vapor  enthalpy  (BTU/lb) 

identity  matrix 

interactive  index 

integral  of  the  absolute  error 

integral  of  the  error 

control  matrices 
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L 

L. 

1 

M 

P 


£ 


s 


S 


T. 

1 

u 


U 

V 

M 

W 

W 

x 


Li 

Vi 

l 


Y. 


i 


z 


— 1  —  3  —  —7 

liquid  flow  from  stage  i  (lb/min) 


-2 

+ 

II  > 

CO 

II  <n 

CO 

<11 

Si 

+ 

ll> 

CO 

G 

s3 

1 — 1 
oil 

+ 

=3 

G 

llo 

CO 

heat  loss  from  stage  i  (BTU/min) 

total  heat  transferred  to  stage  i  (BTU/min) 

reboiler  load  (BTU/min) 

reflux  flow  (lb/min) 

Laplace  variable 

steam  flow  to  reboiler  (lb/min) 

liquid  draw  from  stage  i  (lb/min) 

vapor  draw  from  stage  i  (lb/min) 

temperature  of  stage  i 

T 

control  vector  [QR  D] 

Laplace  transform  of  u 
feedback  control  vector 
discrete  feedback  matrix 
liquid  holdup  on  stage  i  (lb) 
vapor  holdup  on  stage  i  (lb) 
composition  on  stage  i  (wt.  pet.  MEOH) 
vapor  composition  on  stage  i 
(wt.  pet.  MEOH) 

vector  of  states  [x^  , . . .  ,x^q]T 


Z 


Laplace  transform  of  z 
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Greek  Letters 


£ 

fundamental  matrix 

discrete  control  matrix 

A 

discrete  disturbance  matrix 

X.  . 

13 

relative  gain  between  variable  i  and 

variable  j 

Superscripts 


t  _ 

deviation  variable 

steady  state  value 

Of  - 

approximate  case 

A 

discrete  case 

•  — 

time  derivative 

Subscripts 


F 

feed 

FF 

feedforward 

FB 

feedback 

P 

matrix  with  main  diagonal  set  to  zero 

d 

diagonal  matrix 

DC 

decoupling  (for  feedback) 

UC 

decoupling  (for  feedforward) 

' 
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APPENDIX  A 


PHYSICAL  SYSTEM  DATA 

A-l  Schematic  of  Distillation  Column 

Figure  A-l  is  a  detailed  schematic  of  the  pilot  scale 
distillation  column.  In  addition  to  showing  process  streams 
the  diagram  shows  the  DDC  loops  and  multiplexer  points 
associated  with  variables  along  with  the  analog  control 
elements  on  the  column. 

A- 2  Typical  Steady  State  Conditions 

Table  A-l  gives  a  listing  of  the  variables  as  recorded 
during  a  steady  state  run. 

In  Figure  A- 2  some  of  the  major  variables  are  plotted. 

A- 3  Temperature  and  Enthalpy  as  a  Function  of  Concentration 
At  constant  pressure  liquid  enthalpy  is  a  function  of 
concentration.  At  one  atmosphere  the  relation  between 
enthalpy  and  concentration  of  liquid  methanol-water  is  as 
shown  in  Figure  A- 3.  The  relation  used  to  relate  temperature 
and  concentration  is  given  in  Figure  A-4. 

A-4  Methanol-Water  Equilibrium  Curve 

The  equilibrium  curve  for  the  methanol-water  system  is 
shown  in  Figure  A- 5.  Also  shown  is  the  pseudo-equilibrium 
curve  taken  from  experimental  data  obtained  on  the  pilot 
scale  distillation  column. 
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FIGURE  A-l :  A  Detailed  Schematic  of  the  Pilot  Scale  Distillation  Column 
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LPID 

ADC 

DAC" 


HXn 


LEGEND 

indicates  measurements  available  to  IBM  1800 
LPID  -  DDC  loop  identification  code  in  hexidecimal 
-  crosshatching  indicates  no  loop 
ADC  -  multiplexer  input  address  (in  hexidecimal) 
DAC  -  multiplexer  output  address  (in  hexidecimal) 
crosshatching  indicated  no  output 
loop  is  for  data  acquisition  only 

indicates  local  analog  equipment 
TYPE  -  gives  function 
TYPE  CODES 
T  -  transmitter 
R  -  recorder 
C  -  controller 
I  -  indicator 
F  -  flow 

A  -  concentration  (analyzer) 

T  -  temperature 
P  -  pressure 
DP  -  differential  pressure 
n  -  identification  number 

heater  exchanger 
n  -  identification  number 

gas  chromatograph 

thermocouple 

indicates  utilities 
UT  -  ST  -  60  psig  steam 

-  CW  -  cooling  water 

solenoid  valve 
n  -  identification  number 

control  lines 
process  lines 


* 

i  •./  [<  K 
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TABLE  A-l 


STEADY  STATE  DATA 
RUN  NO  SS  3  19/ 


10/73 

FEED  FLOW 

2.400 

LB/MIN 

BOTTOM  PROD 

1.188 

LB/MIN 

REFLUX  FLOW 

1  .936 

LB/MIN 

TOP  PROD 

1.197 

LB/MIN 

STEAM  FLOW 

1.851 

LB/MIN 

COOL  WATER 

37.700 

LB/MIN 

FEED  PLATE 

4 

FEED  COMP 

50.20 

WT  P  C 

TOP  PROD 

97.36 

WT  P  C 

BOTTOMS  COMP 

3.00 

WT  P  C 

FEED  INLET 

161.2 

DEG  F 

REFLUX  INLET 

140.4 

DEG  F 

STEAM  TEMP 

229.7 

DEG  F 

PRESSURE' 

1.4 

IN  H20 

MATERIAL 

B  A  L  A  N  C 

E 

FLOW 

(LB/MIN 

COMP 

)  (WT  PCT ) 

METHANOL 
( LB/MIN) 

WATER 
( LB/MIN) 

FEED  2.400 
BOTTOM  PRODUCT  1.188 
TOP  PRODUCT  1.197 
CLOSURE  ERROR-PC  -0.6 

50.200 

3.000 

97.369 

1.204 

0.035 

1.165 

-0.3 

1.195 

1.152 

0.031 

-0.9 

ENERGY 

B  A  L  A  N  C 

E 

COOLING  WATER 

REFLUX 

TOP  PRODUCT 

FEED 

STEAM 

BOTTOM  PRODUCT 

TOTAL 

HEAT  LOSS 


ENTHALPY  IN 
(BTU/MIN ) 

2298.0 

193.2 

336.4 

2213.4 

5041.1 


ENTHALPY  OUT 
(BTU/MIN) 

3834.7 

198.5 

122.7 

424.9 

248.4 

4829.3 

211.7 


TABLE  A-l  (Continued) 


STEADY  STATE  CONDITIONS  BASED  ON  10  POINTS 
RUN  NO  SS  3  19/  10/73 


FEED  FLOW 

2.400 

REFLUX  FLOW  = 

1.936 

STEAM  FLOW  = 

1.851 

BOTTOM  PROD  = 

1.188 

TOP  PROD 

1.197 

COOL  WATER  = 

37  .684 

TOP  PROD 

97.369 

BOTTOMS  COMP= 

3.000 

FEED  COMP 

50.200 

PRESSURE 

1.458 

COND  LEVEL  = 

5.105 

REB'R  LEVEL  = 

9.953 

DIFF  PRESS  = 

10.887 

REBOILER  TEM= 

208.6 

PLATE  1  TEMP= 

197.3 

PLATE  2  TEMP= 

183.8 

PLATE  3  TEMP= 

174.8 

PLATE  4  T E M P  = 

165.9 

PLATE  5  T EMP  = 

158.2 

PLATE  6  TEMP= 

152.6 

PLATE  7  TEMP= 

149.6 

PLATE  8  TEMP= 

146.7 

COND  TEMP 

143.3 

STEAM  TEMP  = 

229.7 

COND 1 T  TEMP  = 

227.2 

REFLUX  FLOW  = 

124.1 

FEED  FLOW 

92.0 

BOTTOM  FLOW  = 

103.4 

REB  0  ‘  FIE  AD  = 

206.6 

FEED  INLET  = 

161.2 

REFLUX  IMLET= 

140.4 

COL  O' HEAD  = 

147  .4 

WATER  INLET  = 

61.0 

WATER  OUTLET  = 

101 .8 

LB/MIN 

DEV=  0.0061 

LB/MIN 

DEV  =  0.0101 

LB/MIN 

DEV=  0.0100 

LB/MIN 

DEV  =  0.0118 

LB/MIN 

D  EV=  0.0094 

LB/MIN 

DEV  =  0.8833 

WT  P  C 

DEV=  0.0721 

WT  P  C 

DEV  =  0.0000 

WT  P  C  • 

D  EV=  0.0000 

IN  H20 

DEV  =  0.2332 

PSIG 

D  EV=  0.0245 

PSIG 

DEV  =  0.0670 

PSIG 

DEV=  0.2190 

DEG  F 

DEV  =  0.3651 

DEG  F 

DEV=  0.6236 

DEG  F 

DEV  =  0.4533 

DEG  F 

D  EV=  0.5041 

DEG  F 

DEV  =  0.3855 

DEG  F 

DEV=  0.2252 

DEG  F 

DEV  =  0.3405 

DEG  F 

D  EV=  0.2700 

DEG  F 

DEV  =  0.2850 

DEG  F 

D  EV=  0.3259 

DEG  F 

DE V  =  0.2839 

DEG  F 

DEV=  0.3249 

DEG  F 

DEV  =  0.3466 

DEG  F 

D  EV=  0.2312 

DEG  F 

DE  V  =  0.2770 

DEG  F 

D  EV=  0.3313 

DEG  F 

DEV  =  0.2386 

DEG  F 

D  EV=  0.3302 

DEG  F 

DE V  =  0.1900 

DEG  F 

DEV=  0.2124 

DEG  F 

DEV  =  0.5435 

FLOWRATE  (LB/MIN)  COMPOSITION  CHANGE  (WT  PC) 

4*0  -10  1 
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FIGURE  A- 2: 


Experimental  Steady  State  Data. 


LIQUID  TEMPERATURE  (DEG  F) 
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FIGURE  A- 2  Continued 


Experimental  Steady  State  Data 


135-000 


129 


b0 

O 

I — I 

J  o 

<  G 

rG 

Eh  O 
W  <U 
Eh  ^ 
CO  '-n 

W  Mh 

3 

W  cn 

Eh  J  i — I  i — I 

h  Pl,  nj 

Pm  G 


E££"T 


E6T-T 


G 

0) 

+J 

dJ 

3: 

I 


656-0 


o 

G 

fti 

rG 


K 

o 

w 

66Z-0  * 

o 

ss 

Pm 

BBS-0  S 

o 

1-1 

Eh 

< 

EBE-0  h 

S 

P4 

o 

.  2 
o 
o 

66T-0 


000-0- 


Q) 

2 


T3 

0) 

+J 

dJ 

g 

G 

■M 

dJ 

CO 


G 

O 

Mh 


G 

O  • 

•H  0) 
-P  G 
dJ  0) 
G  rG 
+->  ft 
G  co 
<D  O 
O  £ 
G  -P 


0) 

CO  G 
G  O 
CO 

G  -P 
(U  dJ 
> 

CO 

G 

ft  O 
rH  *H 
dJ  -P 
rG  G 
4->  rH 
G  O 
M  CO 


005-0- 


co 

I 

< 

W 

P* 

P 

CD 

M 

Pm 


ENTHALPY  (BTU/LB) 


130 


000 -C 


«PT 


OOVi'OT 


r~' 

b 


H3 

005'/.  § 


000  ’  0 


o 

o 

o 

o 

o 

o 

O 

O 

O 

o 

o 

o 

% 

o 

o 

o 

% 

o 

o 

o 

% 

t 

o 

*r-\ 

ru 

LG 

O 

!TJ 

6 

ru 

31 

V 

m 

a 

u 

0 

•H 

3 

1 

rH 

o 

g 

0 

P 

■M 

0 

£ 

P 

0 

•H 

0 

Ph 

2 
-l-> 
0 
CO 

U 

o 

Pm 


Eh 

G 

O 

o 

P-. 

•H 

+J 

• 

0 

Eh 

Ph 

X 

ooo -a 

V-/ 

■p 

G 

w 

0 

o 

P 

O 

PS 

o 

c 

> 

H 

o 

CO 

Eh 

o 

CO 

V-/ 

< 

W 

PS 

CO  • 

P 

CO 

EH 

3  0 

P 

< 

W 

P 

‘  005“  t? 

P 

M 

CO  Ph 
pH  0 

> 

P 

O 

0  A 

< 

P 

>  Pu 

P 

> 

O 

co 

< 

o 

0  o 

H 

w 

Ph  6 

P 

PS 

0  -H 

W 

p 

-P  < 

£ 

EH 

0 

H 

PS 

<  - 
PS 

■  000 ‘E 

U  0 
0  G 

w 

W 

PUO 

P-l 

Eh 

e 

X 

H 

0  -P 

w 

P 

Eh  0 

4- 

O 

•  • 

H- 

“  005 “T 

1 

< 

W 

PS 

P 

CD 

H 

Pm 


TEMPERATURE  (DEG.  F.) 


WT.  PCT .  MEOH  IN  VAPOR 


131 


FIGURE  A- 5:  Equilibrium  For  Methanol  Water  System  At 

One  Atmosphere. 


APPENDIX  B 


MODEL  EQUATIONS 


B-l  Enthalpy  Model 

Equations  for  non-steady  state  enthalpy  balances  for 
the  individual  stages  of  the  actual  column  are  given  as: 


Reboiler ; 


A-i  =  0 
=  Qr 


h,  =  (h0L0  -  h1L1  -  H1V1  +  Qp  +  q-,)/  CWT1  +  B-,  n  W,n  )  B-l 


ll  '“2  2 
where 


■R 


LI 


11  VI 


VX  =  +  ^1  ”  ^1^ 


L1  =  F  -  D 


Trays  1  to  3;  i  =  2-4 

hi  =  (hi+l  Li+1  +  Hi-1  Vi-1  _  hiLi  _  HiVi  +  QP„ 
(WLi  +  Bli  WVi) 


where 


Vi  =  (Qr  +  ^  qp/ :Hi-hi> 


B-2 

B-3 


B-4 


B-5 


L.  =  F  -  D  +  V. 
l  i-l 


B-6 


Tray  4  (feed  tray);  i  =  5 

h5  =  (hgLg  +  H4V4  +  hFF  -  hgLg  -  HgVg  +  Qg) •  / 
(WL5  +  B15 


(note  that  there  is  only  one  feed  so  the  subscripts  were 
dropped  for  quantities  associated  with  it.) 
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V 


5  =  Qr  +  .Zsl  qj  +  AV 


B-8 


where 


W  A 


Ahp  hp 


and  hp  -  saturated  feed  enthalpy 
Lr  =  F  -  D  +  V„ 


B- 9 


Trays  5  to  8  ;  i  =  6-9 


hi  =  (hi+i  h+i +  Hi-i  h-i  -  hih  -  Hivi  +  ^ 
(WLi  +  BU  wvi) 


B-10 


where 


Vi  =  (Qr  +  F  qj  +  ihFF)/  (Hi  -  hi) 

D  “ 


L.  =  V.  -  D 
l  l-l 


B-ll 


Condenser;  i  =  10 

h10  =  (H9V9  “  h10 

L10 


LI  0 


B-13 

B-14 


B-2  Concentration  Model 

The  material  balance  equations  can  be  written  for  the 
actual  stages  of  the  physical  equipment.  The  expressions 
for  liquid  and  vapor  flow  are  the  same  as  used  in  Appendix 
B-l .  The  resulting  equations  are: 

Reboiler;  i  =  1 

xi  =  Cx2L2  "  X1L1  "  YlV/WLl 


B-15 
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Trays  1  to  3 ;  i  =  2-4 


*i  =  (xi+l  Li+1  +  Yi-1  Vi-1  ”  xi  Li  ‘  Yi  V  /  WLi 


B-16 


Tray  4;  i  =  5  (feed  tray) 


*5  =(X6L6  +  Y4V4  +  XFF  "  X5L5  "  Y5V5)  /  WL5 


B-17 


Trays  5  to  8;  i’ =  6-9 


xi  =  (x.  +  1  l.  +  1  +  Y.^  V._1  -  x±  L.  r  Y.v.)/  WL. 


B-l  8 


Condenser;  i=10 


X10  '  (Y9  V9  ”  x10  (L10  ■  WL10 


B-19 


B-3  Linearization  of  the  Concentration  Model 

In  order  to  obtain  a  linear  model  the  linear  terms  of 
a  Taylor  series  expansion  of  the  concentration  model 
were  taken.  The  general  equation  is 


xi  =  (xi+l  Li+1  +  Yi-1  Vi-1  +  xFi  Fi  Xi  Li 


-  Y.  V.  )  /WT  . 
l  l  /  Li 


B-2  0 


The  linear  version  is  given  by 


x . 

l 


(xi+l  Li+1  +  xi+l  Li+1  +  Yi-1  Vi-1 


+  Yi-1  Yi-1  +  xFi  Fi  +  xFi  Fi  "  xi  Li 
»  »  *  _  /_ 

-  xi  Li  -  Yi  vi  -  Yi  viywLi 


B-2 1 


r'  -  i:?  it*  *  i-iv  x-i Y  4  l*ix)  2  ix 
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The  linearized  vapor  flows  are  given  by 


V.  =  QD  /  (H.  -  h.)  for  i  =  1-4 
1  R/  l  l 

v!  =  (Q*  +  f!  hr .  +  F.  h*.  )/ (H.  -  h.  )  for  i  =  5-10 

l  i\  1  r  l  l  r  l  /l  l 

In  order  to  obtain  a  linear  relationship  between  x  and  Y 
use 


f  i 

Y.  =  C, .  x. 
l  li  i 


B-2  2 


where 


e,  •  = 


6Y . 

l 


'li  <Sx. 

l 


x . 

l 


The  liquid  flow  can  be  obtained  from  mass  balances. 


L.  =  B  +  V.  , 
l  l-l 

for  i  =  1-4 

B-23 

L.  =  B  +  V.  ,-F 
l  l-l 

for  i  =  5-9 

B-24 

L10  =V9  -  D 

B-2  5 

If  this  model  is  put  in  the  state  space  form  given  by 
z  =  £z+gu+Qd  B-2  6 


then  the  elements  of  A  are  given  by: 
Main  diagonal 


a . 
i 


l 


(-E.  -  V.  C,  . 

i  i  li 


w 


Li 


Super  diagonal 


i+1 


B-27 


B-2  8 


■  .  : '  ' 
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Sub-diagonal 


ai,i-l  ‘  (Vi-lCli-lV  ®Li 


B-29 


The  remaining  elements  are  zero 
Letting 


r.  =  (Y.  .  -  x.)  (H.  -  h. )  +  (x.  ,  -  Y. )  (H.  -  h.)  B-30 

1  1-1  1  11  1+1  1  11 


the  elements  of  B  are  given  by: 

for  i  =  1-10 


» .  ..  =  r  .  /  WT  . 
i,l  i/  Li 


B-31 


b .  0  =  (x.  -x.,,)/  W,  .  for  i  =  1-10 
i,2  i  l+l  /  Li 


B-  32 


and  the  elements  of  C  are  given  by: 

ci,i  =  (h+i  "  V/hi  1  =  1-4 


B-  33 


c  •  0  =  c .  0  =  0 

i,2  i,3 


i  =  1-4 


c5 , 1  '  (XF  -x5  +  hF  r5  /  ”L5 


c5 , 2  =  P/WL5 


'/ 


c5 , 3  "  r5  F/WL5 


B-34 
B-3  5 

B-36 
B-3  7 


c  .  n  =  r .  h„  /  WT  . 
i,l  l  F/  Li 


ci  ,  2  =  0 


i  =  6-9 

i  =  6-9 


B-38 


B-39 


c .  0  =  r .  F  /  WT  . 
i,3  l  /  Li 


i  =  6-9 


B-40 


c  =  c  =  c  -  0 

10,1  10,2  10,3 


B-41 
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B-4  Values  for  Concentration  Model 

Table  B-l  contains  numerical  values  for  the  linear 
state  space  form  of  the  concentration  model  using  the 
steady  state  values  as  given  in  Table  B-2.  The  form  of 
the  model  is  as  shown  in  equations  3.10  and  3.11. 

Table  B-3  gives  the  eigenvalues  of  the  model  matrix  A. 


Mj  a  '  <*cl  8«  li»V  S'  ?  54,06  *  lX<fc 


LINEAR  CONCENTRATION  MODEL  MATRICES 
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TABLE  B-2 


STEADY  STATE  DATA 
RUN  NO  • 05D  OL 
8/7/73  0 


FEED 

FLOW 

2.402 

LB/MIN 

BOTTOM  PROD 

1.330 

LB/MIN 

REFLUX  FLOW 

2.215 

LB/M  I N 

TOP  PROD 

1.0  58 

LB/MIN 

STEAM 

FLOW 

1.843 

LB/MIN 

COOL  WATER 

41.897 

LB/MIN 

FEED 

PLATE 

4 

FEED  COMP 

45.60 

WT  P  C 

TOP  PROD 

97.22 

WT  P  C 

BOTTOMS  COMP 

6.18 

WT  P  C 

FEED 

INLET 

162.4 

DEG  F 

REFLUX  INLET 

138.8 

DEG  F 

STEAM 

TEMP 

229.3 

DEG  F 

PRESSURE 

1.3 

IN  H20 

MATER 

I  A  L  B 

ALAN 

C  E 

FLOW 

COMP 

METHANOL 

WATER 

(LB/MIN ) 

(WT  PCT  ) 

(LB/MIN) 

(LB/MIN) 

FEED 

2.402 

45.600 

1.095 

1.307 

BOTTOM  PRODUCT 

1  .330 

6.  180 

0.082 

1.248 

TOP  PRODUCT 

1.058 

97.228 

1.029 

0.029 

CLOSURE  ERROR-PC 

-0.5 

1.4 

-2.2 

ENERGY  BALANCE 


ENTHALPY  IN 

ENTHALPY  1 

(BTU/MIN) 

(BTU/MIN 

COOLING  WATER 

2780.4 

4300.0 

REFLUX 

217.8 

229.9 

TOP  PRODUCT 

109.9 

FEED 

344.4 

STEAM 

2204.4 

422.3 

BOTTOM  PRODUCT 

273.3 

TOTAL 

5547.1 

5335.6 

HEAT  LOSS 

211.4 

TABLE  B-2  (Continued) 


STEADY  STATE  CONDITIONS 
RUN  NO  •  0  5D  OL 


BASED  ON 
8/7/73  0 


FEED  FLOW 
REFLUX  FLOW  = 
STEAM  FLOW  = 
BOTTOM  PROD  = 
TOP  PROD  = 
COOL  WATER  = 
TOP  PROD 
BOTTOMS  COMP= 
FEED  COMP 
PRESSURE 
COND  LEVEL  = 
REB 1 R  LEVEL  = 
DIFF  PRESS  = 
REBOILER  T  EM  = 
PLATE  1  T  EMP  = 
PLATE  2  TEMP= 
PLATE  3  TEMP= 
PLATE  4  TEMP= 
PLATE  5  T EMP  = 
PLATE  6  T EMP  = 
PLATE  7  T EMP  = 
PLATE  8  T E M P  = 
COND  TEMP 
STEAM  TEMP  = 
COND ' T  TEMP  = 
REFLUX  FLOW  = 
FEED  FLOW  = 
BOTTOM  FLOW  = 
REB  O'HEAD  = 
FEED  INLET  = 
REFLUX  I NL  ET  = 
COL  O' HEAD  = 
WATER  INLET  = 
WATER  OUTLET  = 


2.402 

LB/MIN 

2.215 

LB/MIN 

1.843 

LB/MIN 

1.330 

LB/MIN 

1.058 

LB/MIN 

41  .897 

LB/MIN 

97.228 

WT  P  C 

6.180 

WT  P  C 

45.600 

WT  P  C 

1.344 

IN  H20 

5.178 

PSIG 

9.812 

PSIG 

11.734 

PSIG 

206.3 

DEG  F 

188.4 

DEG  F 

175.6 

DEG  F 

169.8 

DEG  F 

166.2 

DEG  F 

157.9 

DEG  F 

152.4 

DEG  F 

149.5 

DEG  F 

147.1 

DEG  F 

144.6 

DEG  F 

229.3 

DEG  F 

226.8 

DEG  F 

126.1 

DEG  F 

96.4 

DEG  F 

110.0 

DEG  F 

204.3 

DEG  F 

162.4 

DEG  F 

138.8 

DEG  F 

147.7 

DEG  F 

66.4 

DEG  F 

102.7 

DEG  F 

10  POINTS 


DEV=  0.0041 
DEV  =  0.0086 
DEV=  0.0146 
DEV  =  0.0083 
DEV=  0.0096 
DEV  =  0.5104 
DEV=  0.0893 
DEV  =  0.0000 
DEV=  0.0000 
DEV  =  0.2553 
DEV=  0.0263 
DEV  =  0.0622 
DEV=  0.4645 
DEV  =  0.3709 
DEV=  0.6800 
DEV  =  0.2838 
DEV=  0.3689 
DEV  =  0.2500 
D  EV=  0.2444 
DEV  =  0.2863 
D  EV=  0.2267 
DEV  =  0.3343 
DEV=  0.2541 
DEV  =  0.2473 
DEV=  0.2610 
DEV  =  0.2582 
DEV=  0.2514 
DEV  =  0.2132 
DEV=  0.1937 
DEV  =  0.3516 
DEV=  0.3837 
DEV  =  0.3395 
D  EV=  0.2037 
DEV  =  0.3456 


■ 

•  ,,  Hi  f  •.  • 
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EIGENVALUES  OF 

1 

2 

3 

4 

5 

6 

7 

8 
9 


TABLE  B- 3 

CONCENTRATION  MODEL  MATRIX 

EIGENVALUE 

-6.495 

-4.500 

-3.523 

-2.769 

-2.278 

-1.150 

-0.973 

-0.652 

-0.244 


10 


-0.055 


APPENDIX  C 


SAMPLE  DERIVATION  OF  CONTROL  LAW  FROM  THIRD  ORDER  MODEL 


In  this  section  a  control  scheme  is  derived  for  a 
general  third  order  linear  state  space  model.  Also  a 
simplified  control  scheme  is  developed  for  comparison. 

The  continuous  state  space  model  can  be  represented 
by 

z  =  Az  +  Bu+Cd  C- 1 

This  can  be  put  in  the  linear  pertubated  form  and  the 
Laplace  transform  taken  to  give 

s-Z  =  AZ  +  BU+CD  C-2 


where 


and 


Z  -  [ z^  Z£  Zg]  -  vector  of  the  states 


1 

U  -  [u-l  UgH  -  vector  of  the  manipulated  inputs 
D  -  d^  or  d  -  disturbance  input 


a.  . 

-  element 

of 

A 

““ 

b.  . 

-  element 

of 

B 

c .  • 

-  element 

of 

C 

— 

then  the  system  represented  by  equation  C-2  can  be 
represented  by  the  signal  flow  graph  in  Figure  C-l. 

If  and  Z2  are  to  be  controlled  using  u-^  and  U2 


respectively,  feedforward  and  decoupling  control  action 
can  be  obtained  by  applying  equation  4.29  to  get 
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s-a  3  3 


FIGURE  C-l :  Signal  Flow  Graph  of  System  Equations. 
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u 


1 


u 


2 


-(al,2Z2 

-(a2,lZl 


+  bl,2U2  +  aI,3Z3  + 

C-3 

+  b2,lUl  +  a2 , 3Z3  + 

C2,ld)/b2,2 

a- 

1 

0 

This  is  shown  by  a  signal  flow  graph  in  Figure  C-2. 

To  obtain  a  solution  these  must  be  solved  simultaneously 
for  u^  and  ^ .  This  is  easily  done  by  analog  elements  but 
for  implementation  on  a  discrete  system  such  as  a  digital 
computer  an  iterative  numerical  solution  is  necessary. 

Another  approach  is  to  use  the  analytic  solution  except 
that  this  does  not  allow  direct  decoupling  of  the  manipul¬ 
ated  variables  as  shown  below.  This  becomes  important 
when  feedback  is  to  be  added  to  the  control  function. 

To  illustrate  this  let  V(z)  be  the  addition  to 
U  to  provide  feedback  action.  Equations  C-3  and  C-4 
become 


U1  =  “ ^ai  2z2  +  bl  2U2  +  al  3Z3  +  C1  ld)//bl,l 

+vl(zl)  C-5 

u2  =  ”^a2  1Z1  +  b2  1U1  +  a2  3Z3  +  C2  l^/^2  2 

tv 2 ( z  2 )  C’6 

where  v^  and  v 2  are  constant  at  any  time  T. 

Solving  equations  C-5  and  C-6  for  u^  and  U2  is 
possible  analytically  but  is  not  practical  because 
the  resulting  function  for  the  manipulated  inputs 
are  functions  of  both  feedback  variables  v^  and 

.  If  these  are  to  be  retained  for  tuning  on-line 
they  cannot  be  incorporated  within  the  function. 
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zi 


Z  2 


Z  3 


pjQURj]  C— 2:  Signal  Flow  Graph  Showing  Control  Relationships 
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The  approximate  solution  from  equation  4.30  yeilds 
the  simplified  feedforward  and  decoupling  control  function 


given  by 


'ui 

1 

II 

b-  1b1 

X  5  X  X  5  C 

0  al,2al,3 

U2 

_  __ 

b2,lb2,2 

a2,l°  a2 , 3 

C-7 


which  has  the  form 


uf 

f  l( Z1  ,Z2  5 z3 

_u2_ 

f  2 ( z-^ , z2  >^2  5 

The  signal  flow  graph  of  this  is  shown  in  Figure  C-3. 
Note  that  the  algebraic  loop  between  u  and  u2  has  been 
eliminated.  The  feedback  V  is  not  decoupled  at  the 
inputs  but  must  rely  on  decoupling  initiated  from  the 


controlled  states. 


/r.iol  Ql  t  8 firf  fta.irfw 
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-  CONTROL  PATH 

-  SYSTEM  PATH 


FIGURE  C-3 :  Signal  Flow  Graph  Representing  Modified 

Control  Equations  With  the  System. 


APPENDIX  D 


CONTROL  MATRICES 

The  control  matrices  used  in  the  simulations  are  given 
in  the  following  figures  along  with  the  discrete  model 
matrices . 

First,  the  tenth  order  continuous  control  matrices  are 
given  in  Figure  D-l. 

Figure  D-2  shows  the  tenth  order  discrete  model.  It 
should  be  noted  that  the  input  matrix  J3  shown  here  contains 
the  coefficients  of  both  control  and  disturbance  inputs. 

The  corresponding  input  vector  is  given  by 

u  =  [F  xp  hp  Qr  D]  D-l 

Figure  D-3  shows  the  tenth  order  control  matrices. 

Figures  D-4  and  D-5  contain  the  fourth  order  reduced 
model  and  control  matrices. 

Figures  D-6  and  D-7  give  the  second  order  reduced 
discrete  model  and  control  matrices. 

Both  the  fourth  and  second  order  model  input  matrices 
contain  coefficients  for  control  and  disturbance  inputs  and 
the  corresponding  input  vector  is  given  by  equation  D-l. 

It  should  be  noted  that  all  the  model  and  control  matrices 
shown  here  are  based  on  the  tenth  order  concentration,  model 
presented  in  Appendix  B. 
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0  2.06E4 

0  1.34E1 


7.16E1 

4.57E-2 


FIGURE  D-l 


MATRIX  Kdc 

0  0  0  0  0  -4 . 4  3E4  0  -6.94E4 
00000  3.36E1  0  5.27E1 

MATRIX  Kff 

0  -1.11E0 

0  8.39E-4 


Continuous  Linear  Tenth  Order  Feedforward  and 
Decoupling  Matrices 


. 


MATRIX 
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FIGURE  D-2:  Tenth  Order  Discrete  Model  Matrices  (T  =  1  Min. 
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FIGURE  D-4:  Fourth  Order  Discrete  Model  Matrices 

(Reduced  from  Tenth  Order  by  Modal  Reduction) 
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FIGURE  D-6:  Second  Order  Discrete  Model  Matrices 

(Reduced  from  Tenth  Order  by  Modal  Reduction) 
(T  =  1  Min. ) 
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FIGURE  D-7:  Second  Order  Discrete  Control  Matrices 

(T  =  1  Min. ) 
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APPENDIX  E 


ON-LINE  COMPUTER  PROGRAMS 

A  series  of  programs  run  on  the  IBM  1800  computer 
enable  control  and  data  acquisition  and  storage  for  testing 
on  the  pilot  scale  distillation  column.  These  programs  were 
originated  by  various  workers  with  some  modifications.  Only 
the  present  usage  of  the  programs  will  be  covered  here. 

A  list  of  the  programs  and  associated  data  files  is 
given  in  Table  E-l. 

In  a  typical  test  the  sequence  of  programs  used  would 
begin  with  the  startup  in  which  DWL  42  would  be  used  to 
activate  the  loops.  Also,  RGM  43  may  be  required  to  reset 
the  shutdown  ECO  which  would  restore  power  to  the  pumps. 

The  startup  of  the  gas  chromatograph  requires  an 
interrupt  on  level  8 ,  bit  6 .  This  is  usually  carried  out 
by  a  hardware  contact.  The  GC  package  queues  a  user 
program  called  DWL  44  which  uses  the  peak  areas  to  calculate 

the  concentration  of  methanol. 

The  program  DASS  is  used  to  in  itiate  a  material  and 
energy  balance  calculation  from  direct  measurements.  It 
links  to  DATAC  which  links  to  BALNC.*  Results  are  used  to 
ascertain  if  the  column  is  at  steady  state. 

Batch  programs  RGM  51  and  RGM  54  are  used  to  enter  and 
write  out  data  contained  in  files  used  by  control  programs. 
These  programs  DWL  52  and  RGM  50  carry  out  the  control  cal¬ 
culation  and  inserts  setpoints  in  DDC  loop  records  as  well  as 
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storing  data  from  DDC  for  up  to  124  minutes  on  a  one  minute 
sample  time. 

It  is  possible  to  commence  another  run  almost  immed¬ 
iately  as  the  data  must  be  copied  to  another  file  using 
RGM  46  for  sorting.  This  frees  the  file  used  to  store  data 

accumulated  during  the  run. 

Finally,  the  data  can  be  printed  and  punched  and/or 

plotted  using  DWL  53. 
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TABLE  E-l 

ON-LINE  PROGRAMS  USED  WITH  DISTILLATION  COLUMN 
DASS,  DATAC,  BALNC  -  do  material  and  energy  balance  from 


direct  measurements  from  column. 

-  use  disk  files  RGM  01  (1  set)  and 

PAC  01  (16  set). 

RGM  4A  - 

writes  labels  to  RGM  01  (for  BALNC  and  DWL  53; 

(batch) . 

DWL  42  - 

turns  column  DDC  loop  on  or  off. 

RGM  43  - 

turns  shutdown  ECO  (re  CRASH)  on  or  off. 

DOUT  (3,16,-) 

DWL  44  - 

queued  by  GC  Pak  to  calculate  MEOH  concentration 

and  put  in  DDC  loop  607.  Uses  CE  Core  words 

Z0017  and  Z001B. 

RGM  46  - 

copies  data  (collected  by  RGM  50)  from  file 

RGM  03  to  RGM  04  for  plotting,  punching,  etc.  by 

DWL  53. 

RGM  51  - 

enters  control  matrices  into  RGM  03  for  RGM  50. 

RGM  54  - 

writes  out  parameters  for  start  run  from  RGM  03. 

DWL  5  2  - 

starts  and  stops  RGM  50  (timer  10)  and  sets 

initial  conditions  for  explorational  run.  Writes 

in  RGM  03  and  to  control  word  in  CE  Core  Z91F9. 

' 
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TABLE  E-l  (Continued) 

RGM  50  -  interrupt  coreload  for  control  and  data 

accumulation.  (Runs  124  min.).  Data  goes  to 
RGM  03  (17  set). 

DWL  53  -  plots,  punches  (prints)  data  transferred  to 

RGM  04  (16  set). 

SSDAT  -  batch  job  to  calculate  internal  liquid  and  vapor 

flowrates  and  vapor  compositions  from  feed,  distri 
bution,  reflux  flows  and  liquid  compositions. 

DWL  50  -  batch  coreload  to  do  thesis  plots. 


APPENDIX  F 
CONTROLLER  CONSTANTS 

The  controller  constants  used  on  the  analog  controllers 
on  the  experimental  equipment  are  given  in  Table  F-l. 

Constants  used  in  controllers  for  the  overhead  composi¬ 
tion  control  tests  in  Chapter  2  are  given  in  Table  F-2. 
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CONTROLLER 

FEED  FLOW 

REFLUX  FLOW 

STEAM  FLOW 

CONDENSER 

LEVEL 

REBOILER 

LEVEL 

PRESSURE 

FEED 

TEMPERATURE 

REFLUX 

TEMPERATURE 

SUPERVISORY 

DDC  LOOP 

FEED  FLOW 

STEAM  FLOW 

DISTILLATE 

FLOW 


TABLE  F-l 

ANALOG  CONTROLLER  CONSTANTS 


PROPORTIONAL 
BAND  (%) 

180 

60 

210 

35 

125 

200 

200 

32 

DDC  LOOPS 

PROPORTIONAL 
CONSTANT  (KP) 

(0601)  -0.9921 

(0603)  -1.4921 

(0605)  -0.4921 


RESET  (MIN.) 
(INTEGRAL) 

0.80 

0.80 

0.15 

8. 

10. 

1. 

0.25 

7. 

INTEGRAL  CONSTANT 
(TI  SEC.) 

89.0 

120.5 

100.0 


I 
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TABLE  F-2 

OVERHEAD  COMPOSITION  CONTROLLER  CONSTANTS 

FROM  CHAPTER  2 


CONSTANTS  FOR  FOXBORO 

CAPACITANCE 

DYNALOG 

PROPORTIONAL 
BAND  (%) 

RESET  (MIN 
(INTEGRAL) 

REFLUX  FLOW  USED  AS 
MANIPULATED  VARIABLE 

200 

7.0 

DISTILLATE  FLOW  USED 

AS  MANIPULATED  VARIABLE 

133 

9.0 

CONSTANTS  FOR  DDC  CONTROL  LOOPS 

PROPORTIONAL 

(KP) 

INTEGRAL 
(TI  SEC. ) 

REFLUX  FLOW  AS 
MANIPULATED  VARIABLE 

1 

o 

• 

CD 

o 

128. 

DISTILLATE  FLOW  AS 
MANIPULATED  VARIABLE 

0.50 

256. 

TUNED 

0.75 

600. 

- 

